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Rozcma,  Daniel  M.,  M.S.,  Purdue  University,  May  1991.  A  Statistical  Analysis  of  a 
Convective  Boundary  Layer  over  Lake  Michigan  oi)  10  January  1984.  Major 
Professor:  Ernest  M.  Agee.  > 

} 

This  research  has  focused  on  the  continued  study  and  analysis  of  data 

I 

collected  over  Lake  Michigan  by  NCAR  research  , 'aircraft  during  a  northerly  How 
cold-air  outbreak  on  10  Jan  84  during  the  Lake  Eufccl  Snow  Studies  (LESS)  projcc.. 
Turbulence  statistics  includi’  g  variances,  covariaiiccs.  TKE,  skewness  and  buoyancy 
have  been  calculated  from  20Hz  data  for  u,v,w,6v  and  q  brr  a  strategically  selected 
30km  flight  segment  and  the  two  contiguous  1 5km  flight  segments  that  comprise  such. 
Emphasis  is  placed  on  the  interpretation  and  comparison  among  the  LESS  results  as 
well  as  those  for  the  Air  Mass  Transformation  Experiment  (AMTEX),  the  Mesoscale 
Air-Sea  Exchange  (MASEX)  Experiment  and  the  models  of  Deardorff  (1980)  and 
Moeng(  19842^  5"^^ 

C_  Of  particular  interest  arc  the  vertical  profiles  ol  normalized  buoyancy  (B) 
and  vertical  velocity  skewness  (Si)  for  the  LESS  results,  where  the  former  shows  a 
primary  maximum  at  level  1  and  a  secondary  maximum  at  level  4  with  the  latter  most 
likely  due  to  condensational  heating.  Jhc  vertical  profile  of  S.v  also  shows  positive 


maxima  at  the  same  levels  indicating  that  skewness  is  buoyancy  driven.  The 
buoyancy  profiles  for  the  AMTEX  and  M.ASFX  botlr  show  only  the  surface 
maximum  with  buoyancy  decreasing  with  height  from  the  surface  to  the  top  of  the 
mixed  layer.  The  buoyancy  profiles  of  c*ase  6  from  the  Deardorff  model  agrees  \  ci  y 
well  with  the  LESS  profile  presented  in  this  sludv. 


1.  INTRODUCTION  AND  STATEMENT  OF  OBJECTIVES 


1 . 1  Literature  Review 

William  Prout  (1834)  is  credited  for  being  the  first  to  introduce  the  term  con¬ 
vection  into  the  field  of  science  as  he  expressed  in  his  own  words;  “There  is  at  present 
no  single  term  in  our  language  employed  to  denote  this  mode  of  propagation  of  heat; 
but  we  venture  to  propose  for  that  purpose  the  term  convection  which  accords  very 
well  with  the  two  other  terms  (conduction  and  radiation)."  The  word  convection 
comes  from  the  Latin  word  convectus  meaning  to  bring  together  or  carry.  Webster 
defines  convection  as  “the  circulatory  motion  that  occurs  in  a  fluid  at  a  nonuniform 
temperature  owing  to  the  variation  of  its  density  and  the  action  of  gravity."  As  a 
result  of  his  pioneering  work  and  thorough  study  while  at  the  University  of  Paris, 
the  name  of  Henri  Benard  (1901)  has  long  been  associated  with  the  study  of  cellular 
convection  that  occurs  in  an  unstable  fluid.  Benard’s  observations  were  anticipated 
to  some  degree  about  twenty  years  earlier  by  James  Thomson  phe  brother  of  Lord 
Kelvin)  who  observed  a  “tessellated  structure”  in  cooling  soapy  water  in  a  tub  seen 
in  the  yard  of  an  inn.  However,  the  most  important  study  in  the  area  of  thermal 
convection  clearly  was  that  of  Benard’s. 

Benard’s  experiments  mark  the  first  quantitative  work  on  the  onset  of  thermal 
instability  and  the  role  played  by  viscosity  in  the  phenomenon.  His  work  consisted 
of  using  a  thin  layer  (approximately  1  mm  thick)  of  whale  oil  called  spermacetti 
that  was  heated  from  below  with  a  horizontal  metallic  plate  maintained  at  uniform 
temperature.  The  fluid  layer,  initially  at  rest,  after  becoming  unstable  resolved  itself 
into  a  number  of  hexagonal  convective  cells  with  ascending  motion  at  the  cell  centers 


and  descending  motion  at  the  periphery  between  adjoining  ceils  (i.e.  uphexagons 
or  closed  cells). 

Lord  Rayleigh  (1916)  went  on  to  theoretically  examine  the  physical  conditions 
leading  to  the  formation  of  convection  cells  in  Benard’s  experimental  study.  Thermal 
convection  (often  referred  to  as  Benard-Rayleigh  convection)  is  recognized  as  the 
systematic  ovenurning  of  a  fluid  layer,  when  sufficiently  and  uniformly  heated  from 
below  and/or  cooled  from  above.  This  overturning  results  in  an  organized  geometric 
pattern  of  convection  cells,  often  viewed  as  he.xagonal  (as  in  Benard’s  study)  although 
a  wide  array  of  geometric  planforms  can  be  achieved  under  a  variety  of  conditions. 
Rayleigh’s  model  for  thermal  convection  represents  a  non-rotating  gravitating  fluid 
with  constant  viscosity  u  and  conductivity  k  initially  at  rest  between  two  free  surfaces, 
and  the  assumptions  of  Boussinesq  (1903)  which  are:  (I)  the  fluctuations  in  density 
which  appear  with  the  ttart  of  motion  result  principally  from  thermal  (as  opposed 
to  pressure)  effects,  and  (2)  in  the  equations  for  the  rate  of  change  of  momentum 
and  mass,  density  variations  may  be  neglected  e.xcept  when  coupled  with  the  gravita¬ 
tional  force.  .Applying  first-order  perturbation  theory  and  the  linearizing  Boussinesq 
approximations  gives  the  following  form  for  the  momentum,  continuity,  state  and 
energy  equations  used  in  the  Rayleigh  model: 
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where  T  is  the  undisturbed  temperature,  a  is  the  coefficient  of  thermal  expansion. 
Pm  is  the  mean  density  of  the  fluid  layer  and  primes  have  been  dropped  from  the 
perturbation  variables  for  convenience.  In  addition,  radiation  and  friction  effects 
have  been  excluded  in  the  energy  equation.  After  nondimensionalizing  the  equa¬ 
tions.  Rayleigh  obtained  three  dimensionless  groupings  of  numbers  that  appeared 
in  the  momentum  and  energy  equations,  each  describing  the  characteristics  of  the 
fluid.  One  of  these  is  the  Rayleigh  number  (Ra,  a  term  coined  by  Chandrasekhar 
(1957))  which  represents  the  ratio  of  buoyancy  forces  and  viscous  forces: 

Ra  =  X5LLi;  (1.7) 

K  V 

Here  g  is  the  acceleration  due  to  gravity,  3  is  the  lapse  rate  of  the  fluid,  d  is  the 
convection  depth,  k  is  the  thermal  diffusivity  and  v  is  the  molecular  kinematic  viscos¬ 
ity.  In  his  theoretical  work.  Rayleigh  showed  that  a  fluid  initially  at  rest  when  uni¬ 
formly  heated  from  below  will  begin  to  convect  above  some  critical  value  of  the 
Rayleigh  number  (Rac)  which  is  dependant  upon  boundary  conditions.  In  Rayleigh's 
model  for  free-free  boundary  conditions,  the  minimum  critical  Rayleigh  number 
(min  Raj  for  which  convection  occurs  is  for  a  Rayleigh  value  of  657.5  at  a  critical 
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wavelength  (kc)  of  2.221. 

Comparisons  between  Lord  Rayleigh's  findings  and  those  of  Benard  revealed 
a  serious  discrepancy  related  to  the  critical  Rayleigh  number  required  for  convection 
to  occur.  The  value  of  the  Rayleigh  number  was  lower  for  convective  onset  in  Be- 
nard’s  experiments  than  it  was  in  Rayleigh’.s  model.  Eventually  Benard  suspected 
that  the  difference  was  perhaps  a  result  of  the  role  played  by  surface  tension  forces 
in  the  convective  process  since  his  experiment  was  with  oil  of  only  a  few  mm  depth. 
Low  and  Brunt  (1925)  recognized  that  the  temperature  gradients  in  Benard's  experi¬ 
ments  were  at  least  tenfold  less  than  that  predicted  by  Rayleigh’s  theory.  Benard 
(1927,  1928)  discussed  the  role  of  surface  tension  as  a  possible  explanation  for 
this  discrepancy. 

As  mentioned  above,  the  critical  Rayleigh  number  for  a  given  fluid  is  depen¬ 
dent  upon  the  boundary  conditions  (free-free,  rigid-free,  rigid-rigid).  Jeffreys 
(1926,  1928)  examined  cases  for  rigid-free  and  rigid-rigid  boundary  conditions 
and  found  successively  larger  values  for  Rac.  For  rigid-free  boundary  conditions 
Jeffreys  found  the  minimum  Rac  to  be  1100.6  and  for  rigid-rigid  conditions  he 
found  it  to  be  1707.7.  These  minimum  critical  Rayleigh  numbers  correspond  to 
critical  wavelengths  of  2.682  and  3.1 17  respectively.  These  results  were  later  con¬ 
firmed  by  Pellew  and  Southwell  (1940)  in  a  linear  convection  model. 
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1.1.1  Geometric  Shape 

Since  a  linear  model,  as  used  by  Rayleigh,  does  not  yield  the  geometry  of 
the  convective  pattern,  Rayleigh  and  his  immediate  successors  were  forced  to  assume 
a  horizontal  planform  for  their  models  that  corresponded  to  an  array  of  square  con¬ 
vection  cells.  By  using  Christopherson’s  (1940)  mathematical  shape  function  for 
a  hexagonal  array,  Pellew  and  Southwell  (1940)  were  able  to  re-examine  Rayleigh’s 
results  for  the  case  of  hexagonal  cells,  as  well  as  for  a  wider  range  of  boundary 
conditions.  The  classical  linear  solution  of  the  relative  vertical  velocity  field  in  the 
hexagonal  cell  obtained  by  Pellew  and  Southwell  is  shown  in  Figure  1.1.  The  zero¬ 
valued  isopielh  is  circular  with  ascending  (or  descending  )  motion  in  the  cell  center 
with  reversed  direction  in  the  periphery  of  the  cell.  Extreme  values  for  the  relative 
vertical  velocities  occur  at  the  corners  of  the  hexagon  and  at  the  cell  center.  Relative 
velocities  at  cell  center  arc  double  those  at  the  corners  of  the  hexagon. 

1.1.2  Transitional  Patterns 

As  pointed  out  by  Agee  et  al  (1973),  the  preferred  convective  mode  is  not 
uniquely  determined  by  the  complete  set  of  equations  of  motion  and  the  boundary- 
conditions  if  the  convening  layer  is  of  infinite  horizontal  extent.  Malkus  and  Veronis 
(1958)  showed  that  for  convective  flows  of  rectangular  and  hexagonal  structure  the 
infinite  number  of  solutions  were  of  a  finite  amplitude.  They  concluded  for  the  case 
of  a  rigid  boundary  and  a  free  surface  that  hexagons  will  not  appear  as  the  initial 
instability.  However,  even  though  vortex  rolls  may  form  initially,  the  tendency  is 
for  a  hexagonal  form  to  evolve  for  a  finite  range  of  disturbances.  Palm  (1960) 
tried  to  explain  this  result  theoretically  and  found  the  variation  in  molecular  viscosity 
with  temperature,  along  with  its  non-linear  effects,  may  be  the  cause.  Segal  and 
Stuart  (1962)  modified  Palm’s  conclusion  saying  that  hexagonal  convection  cells  may 
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be  the  stable  solution  only  when  the  variation  of  molecular  viscosity  with  temperature 
is  sufficiently  large.  As  pointed  out  by  Agee  ct  al  (1973)  two-dimensional  convective 
rolls  (c.g.  cloud  streets)  are  preferred  in  a  convecting  layer  with  significant  vertical 
shear  of  the  horizontal  wind. 

Figure  1.2  shows  a  summary  of  some  of  the  laboratory  and  nonlinear  theoreti¬ 
cal  studies  carried  out  by  Krishnamurti  (1975)  for  a  non-rotating  horizontal  layer 
of  fluid  uniformly  heated  from  below  and  uniformly  cooled  from  above.  The  dia¬ 
gram  shows  the  types  of  flow  observed  and/or  predicted  for  the  ranges  of  the  Rayleigh 
(Ra)  and  Prandtl  (Pr)  number  indicated  where  the  Ra  is  as  defined  above  and  the 
Pr  is  defined  as 

Pr=  V  /  K 

where  v  is  the  kinematic  viscosity  and  k  is  the  thermal  diffusivily.  Along  a  third  coordi¬ 
nate  axis  is  a  schematic  parameter  y  which  is  a  mathematical  measure  of  the  vertical 
a.symmetry  associated  with  any  one  of  a  number  of  physical  processes  that  could  be  im¬ 
posed.  For  Krishnamurti’s  work  the  y  parameter  is  representative  of  the  imposed 
large-scale  vertical  velocity.  The  transition  points  in  the  diagram  depend  upon  factors 
such  as  whether  a  point  is  approached  from  a  smaller  or  larger  Ra  (i.c.  the  hypothesis 
effect) , uniformity  in  the  cell  pattern  (Bussc  and  Whitehead,  1974)  and  imperfections 
such  as  cur\'alurc  and  periodicities  induced  by  the  sidewall  effects  (Krishnamurti, 
1973).  However,  in  general  the  horizontal  convecting  layer  is  seen  to  make  a  slow  tran¬ 
sition  to  turbulence.  At  sufficiently  small  Ra,  the  flow  is  stable  to  all  disturbances.  As 
Ra  is  increased  the  flow  becomes  unstable  to  one  kind  of  disturbance  which  grows  to 
finite  amplitude.  As  Ra  is  increased  even  further,  the  flow  becomes  unstable  to  more 
kinds  of  disturbances.  Once  Ra  is  made  sufficiently  large,  the  flow  is  unstable  to  so 
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Figure  1.2  SJability  regime  diagram  for  Ihc  7  =  0  plane,  where  7  is  a 
measure  of  vertical  asymmedy.  (After  Krishnamurti,  1975) 


many  kinds  of  disturbances  that  the  resulting  flow  is  called  turbulent.  Therefore,  un¬ 
like  the  “fast”  transition  seen  for  Coucitc  flow,  the  horizontal  convecling  layer,  espe¬ 
cially  at  higher  Prantl  number,  remains  in  a  discrete  regime  for  a  finite  range  of  Ray¬ 
leigh  number.  Krishnamurti  (1976)  has  also  derived  the  vertical  velocity  in  hexagonal 
convection,  based  on  nonlinear  theory.  This  hexagonal  solution  is  similar  to  that  by 
Pellew  and  Southwell,  but  it  is  inscribed  within  the  linear  cell  and  rotated  .50'". 

1.1..5  Circulation  Direction 

The  convective  cells  in  Benard's  experiments  had  ascending  motion  at  cell  cen¬ 
ter  surrounded  by  descending  motion  at  the  periphery.  Stommcl  (1947)  postulated  that 
the  favored  direction  of  convective  circulation  involves  the  smallest  frictional  loss,  and 
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that  convective  motion  originates  in  the  low  viscosity  region  toward  the  region  of  high 
viscosity.  Since  the  molecular  viscosity  of  liquids  usually  decreases  with  increasing 
temperature  but  increases  for  gases,  one  would  expect  that  liquids  initially  at  rest 
should  develop  hexagonal  cells  with  ascending  motion  at  cell  center  and  descending 
motion  in  the  periphery  (i.e.  closed  cells)  For  gases,  the  observed  circulation  direction 
should  be  downward  at  cell  center  and  upward  in  the  periphery  (i.e.  open  cells).  It  was 
Hubert  (1966)  who  adopted  the  nomenclature  of  open  and  closed  cells  to  atmospheric 
thermal  convection  ob.scrved  in  convective  marine  boundary  layers.  The  logic  behind 
the  names  comes  from  the  occurrence  of  cloudy  centers  for  closed  cells  and  relatively 
cloud-free  centers  for  open  cells. 

Based  on  the  discussion  given  above,  one  might  expect  to  see  only  open  cells 
during  atmospheric  manifestations  of  Bcnard-Raylcigh  convection.  However  this 
clearly  is  not  the  case.  For  instance,  Rothermel  and  Agee  (1980)  have  documented  the 
coexistence  of  both  open  and  closed  cells  during  the  Air-Mass  Transformation  Ex¬ 
periment  (AMTEX)  over  the  East  China  Sea.  Tippleskirch  ( 19.S6)  attempted  to  inves¬ 
tigate  the  possible  role  played  by  the  vertical  variation  of  molecular  viscosity  on  the  cir¬ 
culation  direction  in  a  Benard-Raylcigh  cell  by  performing  laboratory  experiments 
with  molten  sulfur.  Tippleskirch  cho.se  molten  sulfur  because  it  behaves  like  a  liquid 
in  the  temperature  range  of  1 10°C  to  1-S3®C  with  decreasing  molecular  viscosity  with 
increasing  temperature,  but  behaves  like  a  gas  in  the  range  of  1 5.3 '’C  to  2()0"C  with  in¬ 
creasing  molecular  viscosity  with  increasing  temperature.  In  his  experiment ,  Tipples¬ 
kirch  had  closed  cells  with  the  molten  sulfur  in  the  range  of  1 10'’  C  to  I53'’C  but  ob¬ 
served  a  transition  in  the  circulation  direction  pattern  to  open  convection  cells  in  the 
range  of  153°C  to  200'’ C.  Further  studies  by  Palm  (I960)  and  Segal  and  Stuart  (1962) 
give  additional  support  to  the  idea  that  temperature -produced  changes  in  the  molecu¬ 
lar  viscosity  may  determine  the  direction  of  circulation  in  laboratory  convection. 
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Within  the  atmosphere,  however,  the  transport  of  heat,  momentum  and  mois¬ 
ture  is  dominated  by  the  effects  of  eddy  motions  and  not  molecular  viscosity  and  con¬ 
ductivity.  Based  upon  the  mathematical  solution  of  a  free  convection  model  for  the 
atmosphere  used  in  conjunction  with  Stommcl’s  physical  argument,  Agee  and  Chen 
(1973)  showed  that  the  sign  of  the  change  in  eddy  viscosity  with  height  can  fix  the  direc¬ 
tion  of  circulation  in  cellular  convection.  Or,  more  specifically,  if  the  eddy  viscosity  is 
increasing  (decreasing)  with  height,  an  unstable  layer  of  cellular  convection  will  have 
open  (closed)  cells.  They  go  on  to  emphasize  that  the  mathematical  solution  only 
proves  a  circulation  reversal  if  the  sign  of  the  vertical  gradient  ofeddy  viscosity  changes, 
and  that  Stommel’s  physical  argument  is  necessary  to  determine  the  circulation  direc¬ 
tion.  Krishnamurti  (1975)  later  proposed  that  the  circulation  direction  in  mcsoscalc 
cellular  convection  (MCC)  is  related  to  the  upward  and  downward  large-scale  vertical 
motions  with  closed  (open)  cells  forming  in  regions  of  large-scale  up  motion  (down  mo¬ 
tion).  Van  der  Borght  (1975)  proposed  that  a  stronger  heat  flux  from  the  sea  to  the  at¬ 
mosphere  should  occur  in  regions  of  open  cells  than  in  regions  of  closed  cells.  All  of 
these  studies  suggest  that  many  different  processes  can  affect  or  even  alter  the  direc¬ 
tion  of  circulation  in  a  pattern  of  cells. 

1.1.4  The  Type  I  vs.  Type  11  CTBL 

While  there  are  certainly  more  than  just  two  type ;  of  cloud-topped  boundary 
layers  (CTBLs)  (see  Agee,  1987),  the  Type  I  or  cold-air  outbreak  CTBL  and  the  Type 
11  or  marine  stratocumulus  CTBL  arc  particularly  significant  in  the  earth's  atmo¬ 
sphere.  The  Type  I  CTBL  is  a  thermally-driven  convective  boundary  layer  that  forms 
when  cold  air  moves  over  relatively  warmer  water.  The  Type  1  typically  forms  during 
the  winter  off  the  cast  coast  of  continents  over  warm  ocean  currents  and  south  or  north 
of  the  respective  ice  caps.  Two  particularly  common  regions  for  the  formation  of  the 
Type  I  CTBL  arc  off  the  cast  coast  of  the  United  States  over  the  Gulf  Stream  and  off 


the  East  China  coast  over  the  warm  waters  of  the  Kuroshio  current.  The  two  case  stu¬ 
dies  used  for  comparison  in  this  paper  were  conducted  in  these  areas  and  will  be  dis¬ 
cussed  in  greater  detail  in  Chapter  3.  The  Type  I CTBL  lends  to  form  very  quickly  (  <  24 
hours)  and  may  persist  for  2  -  5  days.  The  cloud  species  formed  in  the  wake  of  the  cold- 
air  outbreak  accompanying  the  Type  1  CTBL  arc  of  the  cumulus  and  stratocumulus 
variety  with  the  geometry  varying  from  2-d  rolls  (i.c.  cloud  streets)  to  2-d  and  3-d 
chains  or  3-d  open  and  closed  MCC. 

The  Type  11  CTBL  is  in  many  ways  the  antithesis  of  (he  Type  1.  It  tends  to  form 
over  water  surfaces  with  temperatures  comparable  to  that  of  thcairabovc.  The  prima¬ 
ry  driving  mechanism  is  cloud-top  radiative  cooling.  Type  II  CTBLs  (end  to  form  in 
the  summer  to  the  west  of  continents  over  cool  ocean  currents  and,  although  they  are 
relatively  slow  to  form,  they  may  persist  for  several  weeks,  with  diurnal  variation  caused 
by  radiative-entrainment  effects.  One  of  the  especially  favored  regions  for  the  forma¬ 
tion  of  the  Type  II  CTBL  is  off  the  west  coast  of  California  as  evidenced  by  expansive 
layers  of  stratocumulus  during  the  summer  in  this  area. 

Both  the  Type  I  and  Type  11  CTBLs  are  characterized  by  weak  wind  shear  and 
negative  large-scale  vertical  velocities  which  act  to  reduce  the  depth  of  the  CTBL.  Dis¬ 
cussion  in  (his  thesis  will  focus  entirely  on  the  Type  1  CTBL. 

Figure  1.3  shows  the  vertical  profile  of  potential  temperature  for  the  typical 
Type  I  CTBL.  The  lowest  layer  of  (he  CTBL  is  the  supcradiabatic  surface  laycrwhich 
results  from  the  very  cold  air  coming  in  contact  with  the  relatively  warm  waters  below. 
Above  this  layer  is  the  well-mixed  laycrwhich  includes  the  neutral  (dry)  subcloud  layer 
and  the  neutral  (moist)  cloud  layer.  Above  this  is  typically  found  the  inversion  layer, 
including  the  transition  layer  through  which  entrainment  occurs,  above  which  is  found 
the  undisturbed  inversion  layer  and  finally  the  free  atmosphere. 


12 


Figure  1.3  Vertical  profile  of  potential  temperature  for  a  Type  1  CTBL  which 
includes  the  transition  layer  where  penetrative  convection  occurs 
(from  Agee  and  Gilbert,  1989). 


1.1.5  Heat  Flux 

Cold  air  outbreaks  (CAOs)  over  regions  favorable  to  MCC  formation  such  as 
the  East  China  Sea  over  the  Kuroshio  Current  and  off  the  East  Coast  of  the  United 
Stales  over  the  warm  waters  of  the  Gulf  Stream  result  in  total  heat  fluxes  (sensible  and 
latent  heat)  that  approach  that  of  the  solar  parameter  (13<St)  W  m  These  large  ener¬ 
gy  fluxes  in  t  urn  arc  capable  of  influencing  large-scale  flows  and  may  provide  a  mecha¬ 
nism  for  providing  energy  for  cyclogencsis  to  occur  (Agee  and  Howlcy,  1977).  Sheuand 
Agee  (1 977)  showed  that  a. sensible  heat  nuxof7t)  Wm  ^and  total  heat  fliixof2tX)Wm  ^ 
were  necessary  conditions  for  the  formation  of  MCC.  They  also  found  that  both  open 
and  closed  cells  occurred  during  the  Air  Ma.s,s  Transformation  Experiment  (AMTEX) 


13 


conducted  over  the  East  China  Sea  when  the  sea  to  air  energy  flux  (latent  plus  sensible) 
was  either  strong  (~1200Wm  -)  or  weak  (~  200  Wm  -).  Previously,  Agee  and  Do¬ 
well  (1974)  had  observed  that  strong  heat  flux  from  the  sea  to  the  atmosphere  occurs 
in  regions  of  open  cells  and  a  weaker  flux  from  the  ocean  to  the  air  in  regions  of  closed 
cells.  Hubert  (1966)  inspected  a  large  number  of  MCC  cases  and  observed  that  open 
cells  can  exist  under  conditions  of  both  strong  and  weak  surface  heating  white  closed 
cells  tend  to  exist  only  with  small  surface  heating.  A  plausible  explanation  for  Hubert’s 
finding  is  that  open  cells  may  be  more  efficient  at  transporting  heat  than  closed  cells 
since  strong  downward  motion  occupies  most  of  the  open  cell  and  continues  to  provide 
cold  air  to  the  sea  surface  and  thus  enhances  the  air-sea  interaction  (Sheu  and  Agee, 
1977).  This  theory  is  supported  at  least  partially  by  the  observance  of  penetrative  MCC 
exclusively  for  the  open  cells;  indicative  of  a  more  vigorous  ascending  branch  in  the 
conveetive  cell. 

1.2  Statement  of  Researeh  Goals  and  Objectives 

Data  collected  over  Lake  Michigan  as  part  of  Project  LESS  have  been  used  in 
this  thesis  research  to  increase  the  unden:tanding  of  convection  in  a  Type  1  CTBL.  In 
particular,  results  from  the  10  January  1984  case  arc  examined  closely  to  study  trans¬ 
port  proce.s.ses  in  a  convective  marine  boundary  layer.  Statistics  including  means,  vari¬ 
ances,  covariances  and  skewness  values  have  been  determined  in  the  hope  of  increas¬ 
ing  the  conceptual  understanding  of  the  structure  of  Type  1  CTBLs.  To  achieve  this 
objective,  the  following  specific  research  goals  were  set  and  have  been  achieved: 

1.  Establish  a  unique  data  set  in  a  homogeneous  region  of  a  Type  I  marine 
CTBL  over  which  turbulent  statistics  can  be  calculated. 

2.  Calculate  the  turbulent  statistics  associated  with  the  heat,  momentum  and 
moisture  budgets  for  a  Type  I  CTBL  (Project  LESS)  using  both  detrended  data  (i.c. 
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w  =  0)  and  raw  data  w  0. 

3.  Compare  and  contrast  the  turbulent  statistics  forihe  It)  Jan  84  ease  of  Project 
LESS  with  other  analogous  obscr\'ational  ease  studies  from  AMI  EX  and  MASEX, 
and  with  model  results  by  Dcardorff  ( 1 980)  and  Moeng  ( 1 984).  with  the  goal  of  explain¬ 
ing  physically  the  similarities  and  differences  of  these  convective  events. 

4.  Investigate  the  role  played  by  the  release  of  the  latent  heat  of  condensation 
in  the  upper  portion  of  the  convective  marine  boundary  layer. 

5.  Compare  .select  conventional  (i.c.  detrended)  data  results  with  non-conven- 
tional  (i.c.  raw)  data  results. 
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2.  PROJECT  LESS 

2.1  Background 

Anyone  who  has  lived  in  Ihe  vicinity  of  the  Great  Lakes  is  well  aware  of  the 
profound  impact  that  these  bodies  of  water  have  upon  the  weather  in  this  region. 
One  of  the  most  obvious  of  these  impacts  is  the  lake-induced  snow  .storms  which 
occur  when  polar  air  masses  move  .southward  out  of  Canada  over  the  relatively  warm 
waters  of  the  Great  Lakes.  This  often  results  in  substantial  increases  in  the  amount 
of  snowfall  at  locations  along  the  downwind  shores. 

Braham  and  Dungey  (198  have  investigated  snowfall  patterns  around  Lake 
Michigan  and  along  the  south  shore  of  Lake  Superior  for  the  winters  of  1909/10 
through  1980/81  with  the  goal  of  quantitatively  estimating  the  effect  of  Lake  Michi¬ 
gan  on  snowfall.  They  have  estimated  the  lake  effect  by  taking  the  ratio  of  observed 
to  predicted  snow  amount  for  four  different  areas  in  Wisconsin  and  Michigan  and 
found  the  long-term  (70  year)  average  effect  along  the  west  shore  (south  of  Sheboyg¬ 
an,  Wisconsin}  was  about  10^8  compared  with  roughly  60%  along  the  east  shore 
(south  of  Hart,  Michigan).  Charnon  (1968)  c.>liiiuiied  lake  effects  on  total  winter 
precipitation  for  the  period  1921  -51  and  found  the  average  effect  along  the  west 
shore  to  be  about  5%  and  and  up  to  30%  in  the  sno\:belt  regions  of  southwest  lower 
Michigan.  Dewey  (1970)  concluded  the  lake  effect  incrca'^es  the  annual  snowfall 
along  the  Wisconsin  shore  by  3(i-40%  while  some  areas  in  western  lower  Michigan 
receive  about  200%  more  snow  than  stations  at  a  comparable  latitude  but  inland, 
aw'ay  from  the  lake  effect. 

The  Lake  Effect  Snow  Studies  (LESS)  Project  was  conducted  during  the  w  inter 
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of  83-84  by  the  University  of  Chicago  Cloud  Physics  Ciroup,  the  National  Center 
for  Atmospheric  Research  (NCAR)  and  the  Purdue  University  Mcsoscnle  Convection 
Research  Group.  The  study  ran  from  1  December  1983  -  28  Januarv  1984.  The 
objectives  of  the  study  included  the  study  of  cold  air  outbreaks  and  resulting  modifica¬ 
tion  of  the  polar  air  that  passed  over  both  Lake  Stiperior  (northerlv  flow  eases)  and 
Lake  Michigan  (both  northerly  and  westerly  flow  cases).  Included  is  an  analysis 
of  convective  PBLs  and  lake-effect  snow  systems.  The  Chicago  Cloud  Physics  Group 
focused  mostly  on  westerly  flow  eases;  the  Purdue  group  on  northerly  flow  eases. 
Routine  ob.scrvations  by  the  Natiorat  Weather  Serv  ice  and  the  Canadian  Atmospheric 
Environment  were  supplemented  by  the  NCAR  Portable  Automated  Mcsonct  (PAM 
11)  and  special  rawinsondc  stations  manned  by  both  the  Chicago  and  Purdue  research 
groups.  Data  were  collected  over  Lake  Michigan  at  roughly  43  N  by  NCAR's  King 
Air  (flight  level  5)  and  Queen  Air  (flight  levels  1  -  4)  aircraft  with  a  total  of  15 
missions  flown.  This  particular  research  has  focused  on  the  data  collected  on  10 
January  1984  between  1812  and  1959  UTC  for  a  northerly  flow  ease.  Figure  2.1 
shows  the  vertical  profile  for  NCAR  aircraft  flights  for  the  10  January  1984  case. 
The  Queen  air  flew  at  flight  levels  1  -  4  while  the  King  -Air  flew  level  5.  From  Figure 
2.1  it  is  noted  that  the  King  Air  flew  through  the  inversion  layer  (marked  by  A), 
the  interfacial  layer  (B)  and  the  mi.xcd  layer  (C)  as  it  flew  from  west  to  east.  The 
Queen  Air  also  flew  throtigh  a  portion  of  the  interfacial  layer  at  level  4.  The  environ¬ 
ment  of  the  interfacial  or  transition  layer  is  characteri/ed  by  warm,  drv  air  from 
the  inversion  layer  mixed  with  the  cool,  moist  air  from  the  convective  region  bv 
entrainment  or  penetrative  convection.  While  these  regions  were  a  focus  of  the  work 
done  by  Agee  and  Gilbert  (1989).  for  the  work  done  in  this  studv  it  was  desirable 
to  obtain  turbulent  statistics  over  as  homogeneous  a  region  as  possible.  This  led  to 
the  selection  of  a  "hand-picked"  data  set  which  not  only  eliminated  as  much  as 
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NCAR  FLIGHT  LEGS  -  10  JAN  84 


20  40  60  80  100  12 


Distance  from  Wisconsin  Shore  (km) 

Figure  2.1  The  five  flight  levels  of  the  NCAR  Queen  Air  and  King  Air 
aircraft  from  1812  to  1959  UTC  for  the  cold  air  outbreak 
of  10  January  1984  (Irom  Agee  and  Mart,  1990) 
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possible  the  effects  of  the  intcrfacial  and  inversion  lavers  hut  also  the  effects  of  the 
Michigan  and  Wisconsin  land  breezes  which  were  detected  at  levels  1  and,  to  a  lesser 
degree,  at  level  2.  This  “hand-picked”  data  set  will  be  discussed  in  more  detail 
in  the  section  on  data  analysis  techniques.  Since  the  10  Jan  84  case  represented 
a  northerly  flow  event,  by  flying  from  west  to  east  (or  vise  versa)  the  aircraft  were 
able  to  make  data  measurements  along  a  path  that  was  transverse  with  respect  to 
the  orientation  of  the  three-dimensional  convective  cloud  bands  present  in  the  flow. 

2.2  Synoptic  Discussion 

CAO  events  are  of  interest  to  the  meteorologist  because  of  the  frequent  dramat¬ 
ic  drop  in  temperature  that  results  from  their  occurrence  and  the  often  limes  severe 
disruptions  to  man’s  activities  in  the  affected  regions.  This  disruption  may  be  espe¬ 
cially  acute  in  the  Great  Lakes  areas  because  of  the  significant  lake-induced  clouds 
and  snow  that  often  times  accompany  such  events.  As  pointed  out  by  Konrad  and 
Collucci  (1989)  such  CAO  events  are  generally  associated  with  a  surface  anticyclone- 
cyclone  couplet  between  which  a  southward  transport  of  cold  air  occurs  from  a  high 
latitude  source  region.  Zishka  and  Smith  (1980)  have  studied  the  climatology  of 
cyclones  and  anticyclones  using  a  2*^  latitude/longitudc  grid  spanning  North  America 
and  surrounding  ocean  waters  for  January  and  July  19.50  -  1977.  They  determined 
areal  distributions  of  cyclones/anticvclonc  events,  genesis,  decay,  rclatixc  variability 
and  preferred  propagation  tracks.  For  the  January  case,  anlicyclogcnesis  generally 
occurs  over  the  snow  and  ice  fickN  from  Montana  through  the  Northwest  Territories 
or  in  the  central  Lfnited  States  from  South  Dakota  to  Texas,  often  south  of  the  jet 
stream.  Using  Wcxlcr’s  classification(1951),  these  two  types  of  anticyclones  appear 
to  be  of  the  “cold”  and  “warm”  anticyclone  class,  respet  tixely.  From  their  analysis, 
Zishka  and  Smith  determined  four  preferred  propagation  tracks  for  January  anticy¬ 
clones  (see  figure  4d  in  their  paper).  One  of  the  preferred  tracks  originates  from 
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the  Northwest  Territories  and  then  splits  into  two  distinct  tracks  over  Manitoba.  One 
branch  continues  southward  and  then  turns  ca.slward  through  the  midwestcrn  Ignited 
States.  The  other  branch  turns  directly  to  the  cast  through  Ontario  and  Quebec. 
The  CAO  event  that  occurred  during  Project  LESS  on  10  Jan  84  was  of  the  first 
type;  a  southward  moving  anticyclone  that  eventually  turned  eastward. 

As  pointed  out  by  Wexler  (1943).  as  a  polar  high  moves  southward,  the  cold 
air  warms  because  of  subsidence  and  surface  heating.  Since  lliis  cold  layer  shrinks 
in  depth  and  covers  an  increasingly  larger  area,  it  typically  decays.  This  was  the 
ca.sc  for  the  CAO  event  of  10  Jan  84.  Behind  the  leading  cold  front  at  12Z  on 
9  Jan  84,  a  strong  high  pressure  (cp  1044mh)  to  the  lee  of  the  Canadian  Rockies 
moved  southeastward  into  Montana  and  North  Dakota.  By  1 2Z  on  10  Jan  84,  the 
cold  anticyclone  was  centered  over  northern  Minnesota  (cp  1 040  mb)  creating  a  light 
to  moderate  north-northwesterly  flow  over  Lake  Michigan.  By  12Z  on  the  11th, 
the  anticyclone  had  moved  southward  over  northern  Illinois  and  had  weakened  an 
additional  8  mb  to  approximately  1032  mb. 

Only  a  brief  account  of  the  synoptic  situation  has  been  given  here;  those 
interested  in  a  more  thorough  synoptic  di.scussion  arc  referred  to  Agee  and  Gilbert 
(1989).  Especially  noteworthy  is  Figure  4  in  the  above  mentioned  article,  which 
shows  a  streamline  analysis  for  both  the  surface  and  850mb  level  flow.  From  this 
figure  it  is  clear  that  the  polar  air  moving  over  Lake  Michigan  analyzed  during  Project 
I.ESS  had  already  been  warmed  as  it  passed  over  the  relatively  warmer  waters  of 
Lake  Superior. 

As  part  of  their  examination  of  extreme  CAOs  over  eastern  North  America. 
Konrad  and  Colucci  (1989)  did  the  following  analysis.  Temperatures  at  the  850mb 
level  from  cold  seasons  (I  December  through  28  February)  during  the  period  1 
January  1976  -  31  December  1985  were  extracted  from  reports  for  20  stations  within 
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the  area  bounded  by  the  30*^  and  50'’N  parallels  and  TS'"  and  90'^W  meridians 
(note  that  the  Project  LESS  operational  area  is  included  in  this  bounded  region). 
The  temperatures  were  then  averaged  over  the  stations  at  1200  UTC  each  winter 
season  day  of  the  10  year  period,  thus  providing  a  regional  mean  temperature  T. 
Results  indicated  a  mean  temperature  for  the  period  in  question  of  -1.7"C  and  a 
standard  deviation  of  5.9'^C.  The  Purdue  Regional  Objective  Analysis  of  the  Mesos- 
cale  (PROAM)  and  its  upper  air  counterpart  (PROAMIJ)  over  this  same  area  gives 
a  rough  estimate  of  the  850mb  temperature  during  the  10  Jan  84  CAO  of  about 
-12°C  (approximately  -15"C  over  the  LESS  area),  nearly  two  standard  deviations 
from  the  mean,  qualifying  this  event  as  a  strong  if  not  extreme  CAO  event. 

2.3  Data  Analysis  Procedures 

As  mentioned  previously  in  section  2.1  above,  for  purposes  of  this  study  it 
was  desirable  to  examine  a  portion  of  the  convectivcly  mixed  layer  over  the  middle 
of  Lake  Michigan  not  affected  by  the  land  breeze  or  inversion  layer  phenomena. 
The  resulting  “hand-picked”  data  set  provides  results  that  can  be  compared  with 
numerical  model  simulations  of  convective  boundary  layers  and  other  relevant  obser¬ 
vational  data  sets  such  as  were  collected  during  the  AMTEX  and  MASEX.  As  noted 
earlier,  the  overall  objective  of  this  study  has  been  to  gain  a  more  thorough  under¬ 
standing  of  the  convective  transport  phenomena  and  associated  turbulence  statistics. 
The  meteorological  parameters  of  interest  for  this  study  include  the  following  20Hz 
data:  u,  v.  w  wind  components  (m/s),  virtual  potential  temperature.  0v  (''K)  and 
specific  humidity  q  (g/kg).  The  data  set  selected  for  this  study  included  data  collected 
from  45-75km  from  the  Wisconsin  shore  at  all  five  verticallv  stacked  levels.  This 
30km  flight  leg  was  in  turn  divided  into  2-151<m  flight  legs  for  which  convergence 
in  the  convective  turbulent  .statistics  was  cheeked.  It  took  1  hour  and  47  minutes 
for  the  NCAR  King  and  Queen  Air  aircraft  to  collect  the  data  used  in  this  study. 
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The  table  below  gives  specifics  about  the  five  flight  legs  flown.  The  heights  indicated 
below  are  meters  above  mean  sea  level  which  arc  easily  converted  to  heights  above 
lake  level  by  subtracting  160  meters.  It  is  clear  from  the  data  shown  below  that 
the  aircraft  did  not  maintain  a  constant  speed  as  they  flew  across  the  lake,  with  flight 
level  5  flown  the  fastest  by  the  King  Air  and  flight  level  1  flown  the  slowest  by  the 
Queen  Air.  Clearly  there  is  more  than  a  sufficient  number  ol  observations  for  each 
flight  segment  to  calculate  meaningful  turbulent  statistics  and  skewness  values. 

Table  2. 1  Characteristics  of  the  flight  paths  flown  by  the  King  Air  and  Queen  Air 
aircraft  for  the  10  January  1984  cold  air  outbreak  during  the  Project  LESS. 


45  -  75  km 

Level 

Height  (m) 

Heading 

5 

1392 

west-cast 

4 

1090 

west-cast 

3 

800 

east-west 

2 

480 

west-east 

1 

210 

east-vtest 

45  -  60  km 

Level 

Heisht  (ml 

Heading 

5 

1392 

west-east 

4 

1090 

west-east 

3 

800 

cast-west 

2 

480 

west-cast 

1 

210 

cast-west 

60  -  75  km 

Level 

HeighLXm) 

Hca.djn_g 

5 

1392 

west-east 

4 

1090 

west-east 

3 

800 

cast- west 

2 

480 

west-cast 

1 

210 

cast-west 

Flight  Legs 


Number 

Speed 

Elapsed 

of_Qiis 

(m/sj_ 

Tirne_Csi 

7969 

75.29 

398.45 

8441 

71.08 

422.05 

8063 

74.41 

403.15 

8256 

72.67 

412.80 

8683 

69.10 

434.15 

Flight  Legs 

Number 

Speed 

Elapsed 

of  Obs 

Lmls} 

Time  (s) 

3938 

76.18 

196.90 

4230 

70.92 

211.50 

4067 

73.76 

203.35 

4143 

72.41 

207.15 

4293 

69.88 

214.65 

Flight  Legs 


Number 

Speed 

Elapsed 

of  Obs 

Time  (s) 

4031 

74.42 

201.55 

421  1 

71.24 

210.55 

3996 

75.08 

199.80 

41  13 

72.94 

205.65 

4390 

68,34 

219.50 
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The  2OH7,  data  for  the  convective  mixed  layer  have  been  detrended  using  a 
second  degree  polynomial  in  order  to  remove  any  synoptic  scale  or  large  mcsoscale 
signal.  Wilczak  and  Businger  (1983)  show  there  is  little  difference  between  linear 
and  quadratic  detrending  when  compared  to  the  variability  when  no  detrending  is 
used.  The  data  were  detrended  using  the  International  Mathematical  and  Statistical 
Library  (IMSL)  subroutine  call  RLFOTH  which  fils  a  univariate  curvilinear  regres¬ 
sion  model  using  orthogonol  polynomials.  The  model  used  for  the  multiple  regression 
is, 

y  =  Bo  +  B|X  B2X^  +  e  (2.1) 

where  y  is  the  value  for  the  given  weather  variable  (u,  v,  w.  0v.  iw  q):  x  is  the  distance 
from  the  Wisconsin  shore  in  meters;  Bo,  B|  and  B2  arc  the  true  coefficients  for 
the  slope  and  x  terms  in  the  equation;  and  e  is  the  random  error,  w'hich  is  assumed 
to  have  a  normal  distribution  with  a  mean  of  7.ero  and  a  variance  of  (t^.  Each  flight 
segment  was  detrended  separately;  i.e.  the  45  -  75  km.  45  -  60  km  and  60  -  75 
km  segments  were  each  detrended  using  a  second  order  model  of  the  form  indicated 
above.  The  value  for  y-hat,  the  estimated  value  for  v  in  the  model  shown  above, 
was  obtained  using  the  IMSL  routine  RLOPDC  which  gives  a  response  prediction 
using  an  orthogonal  polynomial  regression  model.  The  y-hat  value  obtained  for 
each  X  value  was  then  subtracted  from  the  raw  data  values  for  each  variable  to  get 
a  detrended  value  for  that  variable.  The  1  Hz  specific  humidity  data  w'erc  used  to 
calibrate  the  20  Hz  Lyman-Alpha  data  according  to  the  method  of  Friehe  ct  al. 
(1986)  before  they  were  detrended.  The  IMSL  routine  RLDOPM  was  used  to  trans¬ 
form  the  polynomial  regression  model,  fitted  using  orthogonal  polynomials,  into  a 
polynomial  function  of  the  original  independent  variable.  The  resulting  y-intercept 
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and  coefficients  of  the  x  and  terms  in  the  regression  equation  arc  shown  in  Table 
2.2  below. 

Table  2.2  Y-intercepts  and  x-ierm  coefficients  from  2nd  degree  polynomial  de¬ 
trending  process. 

45  -  75  km  Flight  Legs 


Fit  Level 

Intercept 

x-term  coeff 

x^-term  coeff 

u  5 

20.0028 

-7. 08251' -04 

5.2348E-0q 

4 

-9.1543 

3.2419l':-04 

-2.9159E-09 

3 

1.7712 

-7.8785E-05 

8.050SE-IO 

2 

7.5159 

-2.0955I':-04 

1 .5275E  09 

1 

22.7080 

-8.034oF.-04 

6.9054I':-n9 

V  5 

37.1980 

-1.4763F.-03 

1 . 195 5  E-OS 

4 

-20.0785 

4.5535E-04 

-3.4741E-09 

3 

1.6142 

-3.0904E-04 

2.5595I':-09 

2 

-0.8622 

-2.7817E-04 

2.623ir':-09 

1 

-13.1880 

1 .4498E-04 

-S.4877E-in 

w  5 

0.8926 

-3. 1054E-05 

2.3476E-10 

4 

14.3260 

-4.7651E-04 

3.9I39E-09 

3 

5.9528 

-2.  1169K-04 

1 .922  IE-09 

2 

-14.3981 

4.8874E-04 

-4.0054E-09 

1 

-n.S993 

1 .S143E-05 

-1 .4997E-10 

Ov  5 

297. 1423 

-9.5357E-04 

7.2580E-09 

4 

266.0135 

-5.8769E-05 

6.4514F,-I0 

3 

263.2163 

3. 590  IE-05 

-5.2169E-1  I 

2 

261.2700 

9.5538E-05 

-5.25I0E-I0 

1 

263.9518 

8.3481E-06 

1 .4985E-in 

S  ^ 

-8. 1946 

2.7406E-04 

-2.0326E-09 

4 

1.5173 

-1.2223E-05 

1 . 1767E-10 

3 

0.9159 

1.0420E-05 

-6.7952E-1 1 

2 

-0.9655 

8.1903E-05 

-7.0657E-10 

1 

-0.4688 

6,4444E-05 

-5.7129F.-I0 

45 

-  60  km  Flight  Legs 

Fit  Level 

Intercept 

x-term  coeff 

x^-term  coeff 

u  5 

32.2350 

-3.2246E-03 

3.0544E-0S 

4 

71.6853 

-2.7262E-03 

2.5588E-08 

3 

-29.4856 

1 . 1 109E-03 

-1 .0434E-08 

2 

-67.6509 

2.6313E-03 

-2.50O6E-08 

1 

85.7552 

-3.2563E-03 

3.07n3E-ns 

V  5 

43.0201 

-1 .727IE-03 

3. 05 44 E-os 

4 

71.6853 

-2.7262E-03 

2. 5 5831' -08 

3 

-29.4855 

1.  I109E-03 

-1 .0434E.-08 

2 

-6.7650 

2.63I3E-03 

-2.5OO6E-0S 

1 

85.7552 

-3.2563E-03 

3.n7n3r':-ns 

w  5 

12.9967 

-5.  I703E-04 

5.0869E.-09 

4 

54.1998 

-2.0239E-03 

1.SS27E-08 

3 

8.8863 

-3.31  1F.-04 

3.  I325E.-09 

2 

17.4669 

-7. 570  IE-04 

S.0987E.-09 

1 

-3.4638 

I.342ir.-n4 

-1 .3n6H'>nq 

u  5 

20.0028 

-7.0S25E-04 

5. 234SE;-09 

4 

-9. 1543 

3. 24191':- 04 

-2.9I59E-09 

3 

1.7712 

-7.87S5E-05 

S.0508E-10 

2 

7.5159 

-2  0955E-04 

1  .5275E-09 

1 

22.7080 

-8.n34nE-04 

6.9954E,-09 
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Table  2.2,  continued 


Fit  Level 

Intercept 

x-term  coeff 

x^-term  coeff 

Ov  5 

288,7123 

-5.9435K-04 

3,4679E-09 

4 

263,4219 

5.0274K-05 

-2  5892E-I0 

3 

266,7917 

-1.0439E-04 

1  .3l38E-n9 

2 

263.5610 

5.6955E-06 

3.5349E-I0 

1 

254.5918 

3.6696E-04 

-3,2  567E-nQ 

q  5 

-11,7678 

3.9910E-14 

-3.0973E-09 

4 

4.3957 

-1 ,2904E-04 

1 .2950E-09 

3 

3.9764 

-1.  I0I4E-04 

1  . I0S9E-09 

2 

-4.2502 

2.  1418E-04 

-2.0709E-09 

1 

3.3737 

-9.2839E-05 

1  .n666E-09 

60 

-75  km  Flight  Legs 

Fit  Level 

Intercept 

x-term  coeff 

x^-term  coeff 

u  5 

-3.8999 

-1 .2454E-04 

1 .7S37E-09 

4 

103.6111 

-2.91  18E-03 

2.0228E-0S 

3 

89.5051 

-2.6908E-03 

2.ni63E-08 

2 

-82.5159 

2.3678E-03 

-1 .6860E-0S 

1 

28.7860 

-9.9707E-04 

851 15E-09 

V  5 

135.9280 

-4.4188E-03 

3.3784E-0S 

4 

-6.4134 

8.4201E-05 

-1 . 1 121E-09 

3 

15.8231 

-7:2665E-04 

5.6214E-09 

2 

-93.5781 

2.6258E-03 

-1.9881E-08 

1 

-49.0781 

I.2654E-03 

-9.5103E-09 

w  5 

-11.9270 

3.2768E-04 

-2.2679E-09 

4 

29.8591 

-9.3929E-04 

7.3418E-09 

3 

-13.7750 

3.6363E-04 

-2.2592E-09 

2 

-24.5864 

7.7559E-04 

-6.0215E-09 

1 

26. 1639 

-7.6607E-04 

5.5838E-09 

Ov  5 

286.3911 

-5.9921E-04 

4.3983E-09 

4 

277.3638 

-3.9257E-04 

3  2123E-09 

3 

264.7012 

-1.0140E-05 

2.9978E-09 

2 

272.2513 

-2.2991E-04 

1 .8779E-09 

1 

268.3663 

-1 .2622E-04 

1.  17n6E-09 

q  5 

-25.4329 

7.6503E-04 

-5.5169E-09 

4. 

4.9283 

-1 . 1439E-04 

8.7670E-10 

3 

4.5469 

-1 .0228E-04 

7.8075E-in 

2 

1.9937 

-1.3577E-05 

5.5091E-1  1 

1 

18.8223 

-5.  137SE-04 

3. 76  75E-09 

After  the  data  sets  have  been  detrended  for  each  meteorological  variable,  each 
value  represents  a  residual  (i.e.  the  difference  between  the  true  and  predicted  value). 
For  any  regression,  simple  linear  or  multiple,  the  residuals  must  sum  to  zero  regard¬ 
less  of  the  goodness  of  the  fit.  As  a  result,  the  means  for  each  detrended  variable 
must  also  be  zero.  Table  2.3  gives  the  averages  for  all  5  weather  variables  for  the 
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3  different  flight  segments  used  in  this  study.  The  ratio  Z/Zj  is  the  ratio  of  the  actual 
aircraft  height  to  that  of  the  mixed  layer  as  determined  hy  the  mean  height  of  the 
base  of  the  inversion  for  each  flight  segment.  The  mean  inversion  heights  for  the 
45-75  km,  45-60km  and  60-75  km  flight  legs  were  1240m,  1140m  and  1340 
m  respectively.  Note  that  in  the  table,  averages  arc  not  given  at  level  5  for  the 
45-75km  leg  because  this  flight  path  represent  samples  from  both  the  convective 
mixed  layer  and  the  interfacial  layer.  Averages  arc  given  for  the  45-60km  and 
60-75km  legs  because  the  former,  at  least  for  the  most  part,  represents  a  sample 
from  the  interfacial  layer  and  the  latter  from  the  convective  mixed  layer.  Also,  two 
independent  neighboring  samples  allow  one  to  sec  to  what  extent  the  same  results 
can  be  produced  from  within  the  homogeneously  selected  region. 

2.3.1  U  and  V  wind  component  Data 

Figures  2.2  and  2.3  show  the  plots  of  raw  20H7.  U  and  V  wind  components 
(m/s)  for  the  five  vertically  slacked  flight  levels  from  45  -  75  km  with  positive  (nega¬ 
tive)  values  representing  westerly  (easterly)  winds  for  U  and  southerly  (northerly) 
winds  for  V. 

The  plots  of  the  u-component  of  wind  in  figure  2.2  show  that  the  cast-west 
component  was  quite  weak  with  speeds  generally  less  than  3  m/s.  The  strong  land 
breeze  detected  by  Agee  and  Hart  (1990)  extending  from  lower  Michigan  out  to 
around  95  km  at  levels  1  and  2  is  excluded  from  this  “hand-picked"  data  set. 
The  weaker  Wisconsin  land  breeze  that  extended  out  to  about  22  km  from  the  Wis 
consin  shore  at  level  1  is  also  no  longer  evident.  The  plot  of  the  u-component  of 
wind  at  level  5  indicates  a  more  easterly  component  to  the  wind  than  at  any  other 
flight  level.  There  is  also  some  evidence  of  a  transition  from  the  intcrfacial  layer 
(roughly  40-60km  from  the  Wisconsin  shore),  where  the  winds  are  quite  variable. 
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Table  2.3  Weather  variable  averages  for  the  45-75km,  45-  60km  and  60-75  km 
flight  segments. 

Boundary  Layer  Averages 

Averages  (45  -  75  km) 


Flight 


level 

7m 

Avg  1 1 

AvgV 

Avg  W 

Avg  Ov 

Avg  <] 

5 

4 

(Convective  and  non 
0.7500  -0.4188 

-convective  data  not  averaged) 
-6.2302  0.1182 

265.4490 

1.2180 

3 

0.5161 

0.0019 

-7.5212 

0.31.49 

265. 1688 

1.2914 

2 

0.2581 

0.5592 

-7.9129 

0.20.45 

265.0765 

I..4516 

1 

0.0403 

0.2170 

-7.6043 

-0.0220 

265.0095 

1.4172 

Averages  (45  - 

60  km) 

Flight 

l.evel 

7J7a 

Avg  IJ 

Avg  V 

Avg  W 

Avg  0v 

Avg  q 

5* 

1.0807 

-2.3028 

-7.0894 

-0.0.407 

267.1346 

0.6023 

4 

0.8158 

-0.4371 

-6.6398 

0.1900 

265.3432 

1.2149 

3 

0.5614 

-0.1168 

-7..5353 

0.1946 

264.95.38 

1.27.30 

2 

0.2807 

0.9336 

-8.5297 

0.1898 

264,83.48 

1.3806 

1 

0.0439 

-0.0003 

-8.0489 

-0,04(9 

264.8147 

1.4.5.46 

Averages  (60  - 

75  km) 

Flight 

T.evel 

7,/7i 

Avg  U 

Avg  V 

Avg  W 

Avg  0v 

Avgq 

.5 

0.9194 

-4.1473 

-7.8013 

-0.1836 

266.0633 

0.978.4 

4 

0.6940 

-0.4003 

-.5.8188 

0.0461 

265.5553 

1.2210 

3 

0.4776 

-0.1228 

-7..5069 

0.4.453 

265.3876 

1.3102 

2 

0.2388 

0.1820 

-7.2916 

0.2172 

265.3209 

1..4227 

1 

0.0373 

0.4294 

-7.1697 

-0.0022 

265.2000 

1.4816 

Level  5  is  above  Ihe  convective  PBI- 
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Figure  2.2  Plot  of  raw  20Hz  u  wind  component  (m/s)  for  flight  levels 
1  (bottom)  through  5  (top)  with  detrend  curve  overlayed. 
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to  the  convective  mixed  layer,  where,  with  the  exception  of  of  a  segment  from  72 
-  74  km  from  the  Wisconsin  shore,  the  winds  are  steady  at  4  to  5  m/s.  The  u-com 
ponent  of  wind  for  levels  1  through  4  show  much  more  frequent  oscillations  between 
weak  westerly  and  easterly  momentum  than  is  evident  at  level  5.  Levels  1  and  5 
show  the  greatest  synoptic  and  large-mcsoscale  trend.  There  is  also  strong  evidence 
of  a  convergence  band  or  zone  at  approximately  6 1  km  from  the  Wisconsin  shore 
at  level  1  where  there  is  an  abrupt  change  from  westerly  (positive)  to  easterly  winds 
(negative).  It  is  worth  noting  that  the  45-75km  segment  at  flight  level  4.  does  not 
include  the  interfacial  layer. 

The  plots  of  the  raw  v-component  of  wind  in  figure  2.3  show  very  clearly 
that  the  10  Jan  84  CAO  was  a  northerly  flow  event  with  the  northerly  component 
of  wind  typically  2  to  3  times  stronger  than  the  u-component.  The  transition  from 
the  interfacial  layer  to  the  convective  mixed  layer  at  level  5  is  less  evident  in  the 
v-component,  with  the  2-1 5km  segments  appearing  as  mirror  images  of  one  another. 
Interestingly,  the  strongest  northerly  winds  appear  to  occur  over  the  western  portion 
of  the  45-75km  leg  especially  at  flight  levels  1  and  2  (.see  Table  2.3  above).  The 
only  flight  level  for  which  this  docs  not  hold  is  level  5  which  is  influenced  by  the 
interfacial  layer.  Flight  level  4  consistently  had  the  lightest  winds  for  the  three  flight 
segments. 

2.3.2  W  wind  component  Data 

Figure  2.4  shows  the  plots  of  the  raw  20Hz.  w'-componeni  of  wind  for  the 
vertically  stacked  45-75km  segment  flight  paths.  A  very  striking  feature  at  level 
5  in  this  figure  is  the  .strong  peak  in  the  vertical  velocity  at  approximately  57km  from 
the  Wisconsin  shore  marking  the  transition  from  the  interfacial  layer  to  the  convective 


Figure  2.i  Plot  of  raw  lOHz  w  wind  component  (m/V)  for  flight 

levels  I  (bottom)  through  5  (top)  with  detrend  curve  overlayed. 
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mixed  layer.  In  general,  vertical  velocities  at  level  4  are  the  weakest  (1-2  m/s). 
The  average  vertical  velocity  for  level  5  is  negative,  indicating  weak  subsidence  at 
this  level.  For  levels  2,  3  and  4  (within  the  subcloud  and  cloud  layers  respectively) 
the  average  vertical  velocity  is  weakly  positive  (0.2  -  0.4  m/s),  while  for  level  1 

it  is  very  weakly  negative  (-0.002  - 0.04  m/s).  The  accuracy  of  w  and  all  other 

variables  measured  is  presented  in  the  next  section  (see  Table  2.4).  The  w-compon- 
enl  at  flight  level  1  supports  the  possibility  of  a  convergence  band  in  the  vicinity  of 
61km  from  the  Wisconsin  shore  with  strong  upward  motion  (+3.5  m/s)  in  this  area. 
Very  little  trend  is  found  in  the  w-component  of  wind  especially  at  levels  1,  3  and 
5.  However,  each  level  was  still  detrended  using  a  second-degree  polynomial  fit 
for  the  sake  of  consistency. 

2.3.3  Virtual  Potential  Temperature  and  Specific  Humidity  Data 

Figures  2.5  and  2.6  show  the  plots  of  the  20H7,  raw  virtual  potential  tempera¬ 
ture  and  specific  humidity  data  for  the  5  vertically-stacked  levels.  The  intcrfacial 
layer  is  now  clearly  evident  at  flight  level  5  in  both  of  these  profiles. 

For  virtual  potential  temperature  (figure  2.5)  strong  oscillations  (>  3.n'^K) 
are  clearly  evident  in  the  transition  from  the  inversion  laver  to  the  convective  mixed 
layer.  These  oscillations  are  evidence  of  entrainment  spikes  from  the  warmer  inver¬ 
sion  layer  above.  The  data  for  levels  3  and  4  show  very  little  variability  (less  than 
2”K)  across  the  45-75km  flight  path.  Clearly  evident  from  the  overlayed  detrend 
curves  is  the  increase  in  virtual  potential  temperature  as  the  distance  from  the  Wiscon¬ 
sin  shore  increases.  Overall,  the  CTBL  is  weakly  stable  as  indicated  by  the  slight 
increase  in  mean  virtual  potential  temperature  with  height  for  all  flight  segments. 

The  plots  of  raw  20H/,  specific  humidity  data  (figure  2.6)  also  show  the  inter- 
facial  layer  at  level  5  with  very  strong  oscillations  (about  Ig/kg)  over  the  western 


Distance  from  Wisconsin  shore  (m) 


Figure  2.5  Plot  of  raw  20Hz  Virtual  Potential  Temperature  (°K)  for 
flight  levels  1  (bottom)  through  5  (top)  with  detrend  curve 
overlayed. 
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Figure  2.6  Plot  of  raw  20Hz  specific  humidity  (g/kg),  for  flight  levels 
1  (bottom)  through  5  (lop)  with  detrend  curve  overlayed. 
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portion  of  the  flight  leg,  indicative  of  entrainment  of  dry  air  from  the  overlying  inver¬ 
sion  layer.  The  average  specific  humidity  decreases  with  increasing  height  (and  dis¬ 
tance  from  the  moisture  source)  for  all  flight  segments.  The  convective  PBL  appears 
to  be  quite  moist  with  mean  specific  humidities  as  high  as  1.46  g/kg  at  level  1. 
The  specific  humidity  profile  also  supports  the  possibility  of  a  convergence  band  in 
the  vicinity  of  61  km  from  the  Wisconsin  shore  with  relatively  high  (1.6  g/kg)  specific 
humidity  values  in  this  region. 

2.3.4  Detrended  Data 

Figures  2.7  through  2.1  1  show  the  plots  of  the  20H7.  detrended  data  for  u, 
V,  w,  0v  and  q  for  all  five  flight  levels  for  the  45-75km  flight  leg.  A  cursory  compari¬ 
son  between  the  raw  and  detrended  data  plots  shows  that  the  detrending  process  has 
been  successful  in  removing  the  large  scale  trends,  especially  at  level  5,  without  re¬ 
moving  the  microscale  fluctuations.  A  comparison  of  the  maximum  and  minimum 
value  for  each  variable  both  before  and  after  detrending  shows  that  the  latter  has 
very  little  effect  on  the  range  of  the  data. 
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Figure  2.8  Plot  of  detrended  20Hz  v  wind  component  (m/s),  for 
flight  levels  1  (bottom)  through  5  (top). 
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Figure  2.9  Plot  of  detrended  20Hz  w  wind  component  (m/s),  for 
flight  levels  1  (bottom)  through  5  (top). 
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Distance  from  Wisconsin  shore  (m) 

Figure  2.10  Plot  of  detrended  20Hz  vinual  potential  temperature 
(°K),  for  flight  levels  1  (bottom)  through  5  (top). 


Aiiainnii  on.iR  juiainoii  AiKiinnii  :)a.iR 


39 


3.3 

:.s  i 

:.3  • 

■3  i 

3  I- 

-.3  i 
-1.3  i 
-1.3  i 
-3.3  i- 


3.3  r- 

1.3  i 
l.O  i 

.3  I 
0  i- 

-.8  i 

-L.O  '■ 
-1.3  [ 
-3.0  i- 


M  -I  .0 
w  -1.3 
-3.3 


-3.3 

43000 


Distance  from  Wisconsin  shore  (m) 


Figure  2.11  Plot  of  detrended  20Hz  specific  humidity  (g/kg) 
for  flight  levels  1  (bottom)  through  5  (top). 
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2.4  Aircraft  and  Instrument  Specifications 

The  NCAR  Beechcraft  King  Air  (N312D)  and  Queen  Air  (N06D)  are  all 
metal,  low-wing,  pressurized  monoplanes  approved  for  flight  in  known  icing  condi¬ 
tions.  Both  aircraft  have  extensive  research  instrumentation  on  board,  including  sen¬ 
sors  for  measuring  mean  and  fluctuating  components  of  wind  speed  u,  v,  w  and 
u',  v',  w',  ambient  temperature  and  pressure,  dew  point,  water  vapor  fluctuations 
and  geometric  altitude.  Aircraft  position  is  available  from  the  Inertial  Navigation 
System  (INS).  Table  2.2  gives  specifics  about  the  instruments  used  to  collect  weather 
data  used  in  this  study. 

Table  2.4  The  Instrument  specifications  for  the  microphysical  data  collected  by  the 
NCAR  King  and  Queen  Air  aircraft. 


Variable 

Symbol 

Device 

Range 

Accuracy 

Aircraft  Altitude 

ALAT 

INS 

Litton  LTN-Sl 

C- 

1 

<  I.O  nmi 
per  fit.  hr. 

Aircraft  Longitude 

A  I. ON 

INS 

Litton  LTN-51 

-1-/- 1 80'’ 

<  1.0  nmi 
per  fit.  hr. 

Aircraft  Ground 

Speed 

OSF 

INS 

0  -  400  m/s 

<  1.0  knot 
per  fit.  hr. 

Ambient  Temperature 

ATB 

Rosemont 

I02r,2AI. 

-60'’-40'’C 

+/-  .5°C 

Vir.  Pot.  Temperature 

VTIILTA 

Thermoelectric 

I  lygrometer 

-50"  -  50'’C 

-1-/-  L0"C 
<0°C 

Absolute  Humidity 

VLA 

Lyman- Alpha 
Hygrometer 

0,1-  25gtn"^ 

+  /-  5''r  ' 

Specific  Humidity 

SPHUM 

■lohnson-Williams 

LWH 

0  -  5gnr-^ 

Wind  Vector  Hori¬ 
zontal  component 

UI,  VI 

Gust  Probe  and 

INS 

0  -  lOOni/s 
r<10  min) 

+/-  0.  ImA 

Wind  Vector  Vertical 
compotient 

WI 

Gust  Probe  and 

INS 

- 1 5  -  1  5m/s 
(<10  tnin) 

+  /-0.1  tn/s 

'  With  periodic  baselining 
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3.  BOUNDARY  LAYER  STATISTICS 

3.1  Other  Observational  Case  Studies 

As  mentioned  in  Chapter  2,  the  selection  of  a  “hand-picked"  data  set  facili¬ 
tates  a  comparison  of  the  Project  LESS  results  with  those  of  other  observational  case 
studies  and  numerical  models.  Before  the  LESS  results  are  presented,  a  brief  descrip¬ 
tion  is  given  of  these  other  observational  case  studies  (AMTEX  and  MASEX)  and 
the  numerical  models  to  which  comparisons  arc  also  made. 

3.1.1  Air  Mass  Transformation  Experiment  (AMTEX) 

AMTEX  74  and  75  were  conducted  as  a  subprogram  of  the  Global  Atmo¬ 
spheric  Research  Program  (GARP)  sponsored  by  the  Japan  National  Committee  for 
GARP.  The  goal  of  the  field  program  was  to  increase  the  understanding  of  the 
processes  by  which  heat,  moisture,  and  momentum  arc  transferred  from  the  ocean 
surface  through  the  boundary  layer  and  to  the  free  atmosphere  under  conditions  of 
large  temperature  difference  between  the  warm  ocean  and  cold  air  above  (Lcnschow' 
and  Agee,  1976).  The  field  program  was  conducted  during  14-28  February  1974 
and  14  February  -  1  March  1975  over  the  East  China  Sea  over  the  warm  Kuroshio 
Current.  Agee  ct  al.  (1973)  have  indicated  that  this  is  a  climatologically  favored 
region  for  the  formation  of  MCC  of  the  Type  I  variety. 

Sources  of  data  used  during  AMTEX  75  included  33  land  stations  in  Japan 
that  provided  surface  synoptic  data  on  a  regular  basis  during  the  period  of  the  experi¬ 
ment.  Three  special  ships,  Ryqfu  Mont,  Kcifii  Maru  and  Nojima  and  island  "Stations 
were  added  to  the  network  which  formed  a  hexagon  with  Naha.  Okinawa  at  the  center 
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(see  figure  2,  Lenschow  and  Agee.  1976).  Tiie  NCAR  F.lectra  aircraft  flew  nine 
research  missions  originating  fn  n  Naha  totaling  50  hours  4  flight  time.  The  data 
collected  by  the  aircraft  were  the  focus  of  the  analysis  of  Lenschow  et  al.  (1980). 

As  was  done  by  Chou  ct  al.  (1986),  Comparisons  in  this  study  were  restricted 
to  lho.se  between  the  LESS  and  the  AMTEX  results  for  the  four  northcrlv  wind  events 
of  the  AMTEX  that  occurred  on  the  1 5th,  16th.  18th  and  22nd  of  February  1975. 
The  CAO  events  of  the  24th  and  26lh  were  much  weaker,  w  ith  AT  (difference  be¬ 
tween  air  and  sea  surface  temperature)  values  >  -5'^C  (temperature  difference  less 
than  5'^C).  These  AT  values  fall  short  of  one  of  the  necessary  conditions  for  MCC 
formation  proposed  by  Sheu  and  Agee  (1977)  of  a  temperature  difference  of  at  least 
5°C  between  the  overlying  cold  air  and  warmer  water  below.  In  addition  the  fetch 
for  the  4  northerly  wind  events  was  approximately  15b()  km.  w  hile  the  fetch  bu-  the 
northerly  flow  event  of  10  Jan  84  during  the  LESS  was  approximately  250  km. 

.3.1.2  The  Mesoscale  Air-Sca  Exchange  (MASEX)  Experiment 

The  MASEX  field  program  was  conducted  in  January  1983  with  four  daylong 
experimental  flights  conducted  during  CAO  CvCnts  on  16.  18.  19  and  20  Januarv 
1983.  The  primary  goals  of  the  MASEX  experiment  were  to  investigate  the  dynamics 
of  the  marine  PBL  and  to  assess  remote  sensing  techniques  for  measuring  the  ocean 
waves  and  the  surface  sensible  and  latent  heat  fluxes  (Chou  et  al..  1986).  The  area 
of  operation  for  the  experiment  was  bounded  by  35  -41  N  and  76"-69"’\V  as 
shown  in  Figure  3.1  below.  This  area  off  the  east  coast  of  the  Lhiited  States  is  also 
identified  by  Agee  ct  al.  (1973)  as  a  climatologically  favored  region  for  MCC  forma 
tion  as  continental  polar  air  moves  off  the  land  .  ver  the  relativclv  warmer  Gulf  Cur 
rent  waters.  The  experiment  was  conducted  from  the  W  allops  Flight  Facility  with 
three  aircraft,  NOAA  P-3,  NASA  P-3  and  NASA  Electra  collecting  data  via  various 
in  situ  and  remote  sensing  instruments.  Boundary  laver  parameters  were  measured 
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by  the  NOAA  P-3  at  various  flight  levels  within  the  rectangle  ABC'D  shown 


Longitude 


Figure  3.1.  The  NOAA  P-3  and  NASA  Eleetra  aircraft  mcastirement  locations. 
Temperature,  humidity  and  wind  arc  measured  by  dropsondcs  at  locations  S|.  S2 
and  S3.  Downward  looking  lidar  measurements  from  the  N.ASA  Eleetra  aircraft 
arc  marked  by  the  flight  tracks  1  to  5.  Boundary  layer  parameters  are  measured 
by  the  NOAA  P-3  aircraft  at  the  box  ABC  D.  From  Chou  el  al.  fl986) 

in  Figure  3.1.  Temperature,  htimidilv  and  wind  ucre  measured  by  dropsondcs  at 
locations  Si,  S2  and  S3.  Downward-looking  lidar  measurements  from  the  NASA 
Eleetra  aircraft  arc  marked  by  the  flight  tracks  1  to  5.  Fcir  the  northerly  flow  case 
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of  20  January  1983,  the  fetch  was  about  20()km.  Tlie  AT  value  is  not  available 
but  is  probably  larger  than  that  for  either  the  LESS  or  the  AMTEX  because  of  the 
smaller  fetch  and  warm  sea  surface  temperatures  at  the  southeastern  corner 

of  the  box  ABCD).  Table  3.1  gives  a  summary  of  the  characteristics  of  the  two 
observational  case  studies  and  the  resulting  data  sets; 

Table  3.1  Summary  of  observational  case  study  characteristics 


Study 

Fit.  Path 

Elapsed 

Aircraft 

Orientation 

Name 

Dates 

Detrending 

Length  (km) 

Time  (s) 

w'lt  mean  wind 

AMTEX 

15,16,18 

22  Eeb  75 

None 

20  km 

204.8 

Vaiied 

MASEX 

20  Jan  83 

I  .inear 

40  km 

360 

Perpendicular 

LESS 

10  Jan  84 

Quaciratic 

30  km 

Varied 

300-430 

Perpendicular 

The  dates  indicated  for  the  AMTEX  are  not  inclusive  of  all  the  dates  of  the  study 
but  indicates  only  those  used  for  comparison  in  this  paper  (where  air-sea  temperature 
differences  exceeded  5^C). 

3.2  Numerical  Models 

Turbulent  statistics  calculated  frorr  the  Project  LESS  data  were  also  compared 
with  those  from  (wo  models;  Deardorffs  3-D  model  (1980)  and  Moeng’s  Largc-Ed- 
dy-Simulation  model  (1984).  Each  model  is  briefly  discussed  below. 

3.2. 1  Dcardorff’s  3-D  Model 

The  Deardorff  model  was  a  three-dimensional  numerical  model  containing 
40’  grid  points  inside  a  volume  of  horizontal  area  5km  x  5km  and  height  2km. 
The  horizontal  and  vertical  grid  spacing  were  125m  and  50m  respectivelv.  The 
situation  treated  was  day-33  of  the  Wangara  data  between  9a.m.  and  5p.m.  The 
surface  was  assumed  to  be  horiz.ontallv  homogeneous  with  respect  to  temperature 
and  roughness.  The  scaling  parameter  Zj,  the  depth  of  (he  mixed  layer,  was  defined 
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by  Deardorff  as  the  height  at  which  the  sensible  heat  flux  is  most  negative.  For 
purposes  of  comparison  only  case  6  from  Dcarclorff’s  work  will  be  used  since  it  most 
closely  replicates  the  conditions  of  the  LESS.  This  case  is  for  a  stratocumulus  deck 
occupying  the  upper  half  of  the  mixed  layer  with  liquid-water  effects  included  and 
cloud-top  radiative  flux  divergence  present  at  the  uppermost  grid  volume  in  each 
vertical  column  for  which  qi,  the  Reynolds  average  of  specific  liquid-water  content 
over  the  volume  Ax  •  Ay  •  Az,  exceeded  0.01x10'-  (Deardorff.  1980).  For  this 
case  the  physical  domain  was  reduced  to  a  2km  cube  and  the  horizontal  grid  spacing 
was  set  at  50m.  The  radiative  flux  divergence  then  occurred  over  only  the  uppermost 
50m  of  cloud,  while  the  cloud-top  height  varied  somew'hat  in  space  and  time. 

.^.2.2  IVfocng’s  Largc-Eddy-Simulalion  Model 

As  pointed  out  by  Moeng  (1984),  a  model  which  explicitly  simulates  the  large 
eddies  but  parameterizes  the  small  eddies  is  called  a  large-eddy  simulation  (LES) 
model.  Moeng’s  LES  model  used  a  pseudospectral  representation  in  (he  horizontal 
direction,  in  which  the  PEL  turbulence  is  nearly  homogeneous,  and  second-order 
finite-differencing  in  the  vertical  direction.  The  model  output  (to  be  discussed  later) 
showed  that  the  resolvable-scale  eddies  contain  more  turbulent  energy  and  transfer 
more  heat  than  do  the  subgrid-scalc  eddies,.  The  latter  act  mostly  to  dissipate  the 
turbulent  energy  which  cascades  down  from  (he  resolvable-scale  eddies.  Except  for 
within  the  surface  layer,  (he  subgrid-scalc  eddies  do  not  contribute  much  to  direct 
or  indirect  heat  transfer  (Moeng,  1984).  The  l.ES  model  was  used  to  simulate 
day- ."^3  of  the  Wangara  data  from  hours  9  and  16.  The  grid  used  measured  32 
X  32  X  40  resulting  in  a  horizontal  and  vertical  resolution  of  approximately  156m 
and  50m  respectively  with  a  3s  time  step  used.  For  simplicity,  the  moisture  field 
and  radiation  effects  were  left  out  of  the  model. 
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3.3  Scaling 

In  this  and  the  following  sections,  boundary  layer  statistics  will  be  examined 
for  the  10  Jan  84  CAO  event  of  the  Project  LESS.  Included  in  the  statistics  are 
means,  variances,  covariances,  various  TKE  component  values  and  skewness  values. 
Results  are  presented  for  the  45-75km  flight  segment  and  the  2-1 5km  (45-60km 
and  60-75km)  segments  which  comprise  it.  In  accordance  with  similarity  principles, 
appropriate  scaling  has  been  done  for  the  turbulent  statistics  based  on  Stull  (1988) 
and  Lenschow  et  al.  (1980)  as  follows: 

Zj.'  depth  of  the  mixed  layer  above  the  ocean  surface,  taken  to  he  the  inversion 
base  height; 

w*:  free  convective  scaling  velocity,  based  on  the  buoyancy  flux  at  flight  level  1 
(50m  above  the  lake  surface); 


w 


*  _ 


(w'©'v)i 


u*:  friction  velocity,  U* 


[ 


-  2  -  2 

(  u'w'  )|  +  (v'w'  )i 


] 


(3.1) 


(3.2) 


(]■'':  moisture  scale. 


q*  =  (w'q'), 

w  * 


(3.3) 


©v  :  virtual  potential  temperature  scale,  0^*  =  (  \v'©v')| 


(3.4) 


w 
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The  free  convective  scaling  velocity  has  been  used  for  scaling  in  most  instances  be 
cause  the  boundary  layer  was  well  mixed  and  filled  with  thermals  due  to  strong  heat¬ 
ing  from  below.  The  friction  velocity  has  been  used  for  scaling  the  stress  profiles 
so  as  to  facilitate  comparison  with  the  MASEX  and  AMTEX  results. 

3.4  Boundary  Layer  Means 

Figure  3.2  gives  the  vertical  profile  of  means  for  the  raw  20Hz  u  and  v  wind 
components  for  the  45-75km,  45-60km  and  60-75  km  flight  segments.  Note  that 
mean  wind  speeds  are  not  indicated  for  flight  level  5  for  the  45-75 km  leg  since 
this  level  represents  a  mixture  of  both  convective  and  non-convective  data.  Level 
5  for  the  45-60  km  leg  is  within  the  interfacial  layer.  This  is  particularly  evident 
in  the  u-component  which  has  its  highest  mean  value  at  this  level.  The  overall  mean 
wind  is  almost  due  north  throughout  the  mixed  layer  with  speeds  of  between  1 2  and 
15  knots.  For  the  45-60  km  segment  at  level  5  the  mean  wind  is  at  20”  (i.c.  from 
the  northeast)  at  about  55  ms"*  (approximately  15  knots).  The  strongest  average 
winds  are  in  the  subcloud  layer  at  levels  1  and  2;  the  weakest  are  in  the  cloud  layer 
at  levels  3  and  4. 

Figure  3.3  (see  below)  shows  the  vertical  profile  of  the  means  for  the  raw 
20F17.  vertical  velocity,  w,  for  the  different  flight  segments.  Clearly  the  mean  vertical 
velocity,  w,  is  not  zero  at  any  of  the  flight  levels  as  is  required  in  many  numerical 
models.  The  strongest  (positive)  mean  vertical  velocities  arc  at  level  3,  within  the 
cloud  layer,  for  all  three  segments;  the  weakest  arc  at  level  1  where  the  mean  vertical 
velocities  are  weakly  negative.  The  mean  vertical  velocity  at  level  5  for  both  15km 
segments  is  weakly  negative,  indicative  of  overall  weak  'subsidence  at  this  level. 

Figure  3.4  (sec  below)  is  the  vertical  profile  of  the  means  for  both  the  raw 
20Hz  virtual  potential  temperature,  0v,  and  the  raw  2()Hz  calibrated  specific  humid¬ 
ity,  q.  As  one  would  expect  in  a  Type  F  convective  boundary  layer,  the  0v  profile 
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Figure  3.2  Vcrlical  profile  of  mean  raw  2()f^/  ii  and  v-wind  speeds  (m/s) 
for  the  5  ProjccI  LESS  flight  levels  for  the  45-75km, 
45-r>nkm  and  60-75km  segments 


Figure  3.3  Vertical  profile  of  average  raw  20Hz  w-wind  speeds 
for  the  5  Project  LESS  flight  levels  for  the 
45-75km,  45-60km  and  60-75km  segments 
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is  neutral  dry  within  the  subcloud  layer  (level  2)  and  neutral  moist  within  the  cloud 
layer  (levels  3  and  4).  The  vertical  profile  of  0v  shows  very  clearly  the  warm  air 
intrusion  due  to  entrainment  occurring  at  the  top  of  the  PBL  for  both  the  45-60km 
and  60-75km  segments  with  an  increase  in  0v  of  nearly  2"K  from  level  4  to  level 
5  for  the  45-60km  flight  segment.  The  vertical  profile  of  q  shows  a  decrease  in 
specific  humidity  w'ith  height  as  the  distance  from  the  moisture  source  (i.e.  Lake 
Michigan)  increases.  The  drop  in  q  is  especially  dramatic  for  the  45-60km  segment 
at  level  5  where  a  50%  decrease  is  observed  from  level  4  indicating  entrainment 
of  dry  inversion  air. 
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3.5  Boundary  Layer  Variances 

As  pointed  out  by  Stull  (1988)  variances  of  variables  such  as  u,  v.  w,  0v, 
and  q  give  information  about  the  turbulent  energies  and  intensities.  Lcnschow  and 
Wyngaard  (1980)  point  out  that  turbulence  second-moment  budgets  are  of  consider¬ 
able  importance  in  atmospheric  boundary  layer  (ABL)  modeling.  In  this  section 
the  average  variances  of  u,  v,  w,  0v,  and  q  are  prc.sentcd  as  computed  from  the 
detrended  data  sets  obtained  from  the  raw  20H7.  data  as  described  in  section  2.3. 
To  calculate  the  turbulent  variation  for  a  variable,  u,  one  would  normally  subtract 
the  mean  u  from  each  instantaneous  value  of  u  as  shown  in  the  equation  below; 

u'  =  Uj  -  u  ,  (3-5) 

where  u'  represents  the  turbulent  fluctuation  in  u.  This  process  is  known  as  Reynold's 
decomposition.  For  the  detrending  process,  u  can  he  thought  of  as  a  running  average 
(actually  its  the  predicted  value  for  each  Uj  resulting  from  the  second  degree  polyno¬ 
mial  detrending)  which  is  subtracted  from  each  instantaneous  value  of  u.  The  result¬ 
ing  u'  values  are  squared  and  then  averaged  over  the  given  fMght  segment.  Finally, 
the  averages  are  normalized  by  the  square  of  the  appropriate  scaling  factor  (w*,  0v* 
or  q*)  which  provides  dimensionle.ss  results  for  comparison  with  other  studies. 

Table  3. 1  gives  a  summary  of  (he  variances  for  the  3  flight  segments  for  the 
10  Jan  84  CAO  event.  The  values  for  w*,  0v*,  and  q*  are  also  shown  for  each 
flight  segment.  Since  each  segment  was  detrended  separately,  the  scaling  factors  were 
also  calculated  separately  for  each  leg. 
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Table  3.2  Project  LESS  boundary  layer  normalized  average  variances. 


Variances  (45  -  75  km) 


Zi  =  124()m,  w*=  1.4684  ms  0,.*  =  0.0470  ” K,  q*  =  0.02SI  g  kg  ' 


Flight 

Level 

Z 

Zi 

W^2 

XV  *2 

Zll 

©Z! 

C 

5 

(Convcclivc  and  non-convcclivc  data  not 

averaged) 

1.8.536 

4 

.7500 

0.5097 

0.4377 

0.5700 

5.7031 

3 

.5161 

0.4267 

0..3666 

0.5669 

2.0492 

7.0896 

2 

.2581 

0.5047 

0.71.32 

0.6588 

.3.5537 

10.7675 

1 

.0403 

0..3714 

0.3587 

0..5452 

6.4843 

I.5.2I17 

Variances  (45  -  60  km) 


Zi  = 

1140m,  w*  = 

1.4445  ms 

‘,0V*  =  0.0495  '^K, 

q*  =  0.0.347  g 

!  kg 

Flight 

Z 

ELI 

Level 

Zi 

W*2 

W*2 

0V 

9 

5* 

1.0807 

0.5294 

0.4644 

0.2090 

185.8579 

37.2.395 

4 

.81.58 

0.5802 

0.4676 

0.6.392 

3.4100 

2.2073 

3 

..5614 

0.47,53 

0.3125 

0.5701 

1.8690 

4..3771 

2 

.2807 

0.4269 

0.5486 

0.6)92 

3.4206 

,3.7420 

1 

.04.39 

0.3252 

0.3.394 

0..5907 

5.7819 

5.8629 

Variances  (60  -  75  km) 


Zi  = 

1340m,  w'*  = 

1 .4754  ms 

1, 0.0440  "K, 

q'^  =  0.0194  g  kg  ' 

Flight 

Level 

Z 

Zi 

W*2 

vviz 

©v'  ^ 
0v*2 

5 

.9194 

0.6318 

0.4262 

0. 1 .390 

37.0043 

17.6242 

4 

.6940 

0.2210 

0.,3606 

0.4828 

2.9,577 

1.8473 

3 

.4776 

0.3227 

0.3840 

0.5723 

1.8727 

7.4173 

2 

.2.388 

0.3299 

0.3478 

0.683 1 

2.7842 

29.3327 

1 

.0.373 

0..3.162 

0.2983 

0.5085 

5.5732 

15.7693 

^  I  x’vcl  5  is  above  I  lie  convective  I’ltl , 
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Note  that  the  free  convective  scaling  velocity  is  generally  4  to  5  times  as  large  as 
the  friction  velocity  (see  Table  3.3  for  u*),  indicating  that  buoyancy  is  dominating 
over  shearing  forces  at  level  1 . 

3.5.1  Comparisons  with  Observational  Case  Studies 

Figure  3.5  shows  the  vertical  profile  of  the  average  u  variance  for  both  the 
Project  LESS  and  AMTEX  studies.  Generally  the  shapes  of  the  profiles  are  in  very 
good  agreement  with  values  for  the  LESS  approximately  500'  to  100%  larger  than 
the  AMTEX  values  with  the  exception  of  the  very  lowest  data  levels  where  the  AM¬ 
TEX  values  are  slightly  larger.  The  AMTEX  results  may  be  higher  than  for  the  LESS 
at  the  lowest  levels  because  data  were  collected  closer  to  the  ocean  surface  where 
shear  would  be  strongest.  Since  the  mean  variances  have  been  scaled  using  vv*. 
which  takes  into  account  the  surface  buoyancy  flux  and,  hence,  the  different  strengths 
of  the  heat  source  (i.e.  different  AT  values),  the  larger  values  of  u  variance  for  the 
LESS  is  probably  due  to  .stronger  shear  resulting  from  a  shorter  fetch  than  that  of 
the  AMTEX  (250km  vs  1500km) 

The  vertical  profile  for  the  average  v-component  variance  (figure  3.6)  is  very 
similar  to  that  of  the  u-component.  The  values  for  the  LESS  are  consistently  larger 
than  those  of  the  AMTEX;  as  much  as  300%  larger  for  level  2  of  the  LESS  data. 
Very  little  variation  is  noticeable  in  the  AMTEX  results;  but  a  significant  peak  is 
clearly  shown  at  level  2  in  the  LESS  profile,  with  a  value  of  0.71  32  for  the  45-75km 
leg.  The  vertical  profile  for  the  average  w-component  variance  (figure  3.7)  shows 
more  clearly  than  either  of  the  first  two  profiles  how  much  more  dynamically  active 
the  CTBL  for  the  LESS  was  than  that  of  the  AMTEX.  Values  arc  consistently  200% 
to  300%  larger  in  the  LESS  case  with  maxima  found  in  the  subcloud  layer  at  level 


.5 
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Figure  3.5  Vcriical  profile  of  average  u-wincl 
componenl  variance  for  the  LESS 
and  AMTEX  normalized  hv  u 
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Figure  3.6  Vertical  profile  of  average  v-wincl  component  variance 
for  the  LESS  and  the  AMIEX  mninalized  hy  vv*^ 
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Figure  3.7  Verlical  profile  of  average  w-\vind  coinponenl  variance  for 
the  LESS  and  (he  AMTEX  normalized  by 
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2  and  in  the  cloud  layer  at  level  4.  Both  the  LESS  and  AMTEX  results  show  the 
effect  of  the  capping  inversion  with  a  significant  drop-off,  especially  in  the  LESS 
results,  in  the  average  variance  at  levels  near  the  inversion  base.  Clearly  the  capping 
inversion  acts  to  dampen  the  vertical  velocity  variances  (i.c.  turbulence)  at  this  level. 

The  vertical  profiles  for  the  average  0v  and  q  variances  (see  Eigures  3.8 
through  3.1  1)  show  the  effects  of  entrainment  verv  clearly  near  the  top  of  the  mi.xcd 
layer,  where  the  variances  arc  as  much  as  50  times  larger  than  over  the  bulk  of  the 
mixed  layer.  The  increase  in  the  0^.  variance  near  the  top  of  the  mixed  layer  is 
larger  than  that  for  q  variance  because  the  normalized  potential  temperature  jump 
(A0v/0v*)  is  larger  than  that  of  the  humidity  jump  (Aq/q’'  )  at  the  top  of  the  PBL 
(Chou  ct  al.,  1986).  Results  for  the  MASEX  arc  from  segments  of  the  box  ABCD 
shown  in  Figure  3.1  earlier.  Both  the  AB  and  CD  flight  segments,  both  of  which 
were  perpendicular  to  the  mean  wind  direction,  are  divided  into  an  cast  (E)  and 
west  (W)  segment.  This  division  was  done  by  Chou  ct  al.  because  the  PBL  height 
increases  eastward  and,  hence,  the  normalized  height  Z/Z|  decreases  eastward  for 
a  leg  flown  at  a  constant  height. 

3. .5. 2  Comparisons  with  Model  results 

A  comparison  of  (he  average  u  variance  for  the  1-ESS  w'ith  results  from 
Moeng’s  LES  model  for  hour  1  6  of  day  33  of  the  Wangara  experiment  (Figure  3. 1  2) 
reveals  a  very  similar  profile  in  the  lower  half  of  the  convective  PBL.  In  the  upper 
half  of  the  PBL  the  u  variance  for  the  l-ESS  continues  to  increase,  reaching  a  maxi¬ 
mum  at  level  5.  In  the  LES.  the  average  u  variance  decreases  above  0.5  Z/Z|. 
The  increase  in  average  u  variance  for  the  LESS  may  be  due  to  the  release  of  the 
latent  heat  of  condensation  at  level  4,  thus  increasing  (he  horizontal  component  of 


TKE. 


Figure  3.10  Vcrtieal  profile  of  coverage  \irlual  poteruial 
icmperalure  variance  for  the  f.ESS  and  ihc 
MASF.X  normalized  bv 


Figure  3.1  1  Vertical  profile  of  average  specific  humidity  variance 
for  the  LESS  and  the  MASEX  normalized  hv 
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3.6  Boundary  Layer  Covariances 

As  pointed  out  by  Stuil  (1988),  covariances  piovidc  a  measure  of  flux  or 
stress,  where  a  flux  is  defined  as  the  transfer  of  a  quaniiiy  such  as  mass,  heat, 
moisture  or  momentum  per  unit  area  per  unit  time.  Statistically,  the  covariance  of 
two  variables  A  and  B  is  defined  as; 


N-1 


■Z 


covar(A,B)  =  1  /  (Ai-A)‘(B,-B) 

N 

i=0 


(3.6) 


Using  Reynolds  averaging,  this  can  be  reduced  to; 


N-1 


covar(A,B)  =  1 


N  i=0 
=  a'  b' 


a'i  b'i 


(3.7) 


The  covariance  then  indicates  the  degree  of  common  relationship  between  the  two 
variables  A  and  B.  The  covariances  that  have  been  calculated  and  will  be  presented 
in  this  paper  arc  those  listed  below: 

Average  vertical  kinematic  eddy  flux  of  u-momentum  =  u'w' 

Average  vertical  kinematic  eddy  flux  of  v-momenium  =  v'w  ' 

Average  vertical  kinematic  eddy  flux  of  virtual  potential  temperature  =  u'0'v 
Average  vertical  kinematic  eddy  moisture  flux  =  w'q' 

Future  references  to  these  quantities  wall  omit  the  w'ord  kinematic. 

3.6.1  Comparisons  among  Observational  Case  Studies 

Table  3.3  shown  below  gives  fluxcs/slrcsses  for  the  buoyancy,  moisture  and 
momentum  for  the  Project  I. ESS  data  set.  As  was  done  previously  for  the  momcn 
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lum,hcat  and  moisture  variances,  the  covariances  have  liecn  normalized  by  the  qiian 


lilies  indicated. 

Table  3.3  Project  LESS  boundary  layer  normalized  average  covariances 

Covariances  (45  -  75  km) 


Zi  =  1240m,  vv*  = 

1.4684  ms  ' 

0.3159  ms  '. 

,  0/ =0.0470 

K.  q' =  0.0281  gkg 

Flight 

z 

w'q' 

vv'0v' 

uA\ ' 

a. 

V  w' 

Level 

Zi 

w*q* 

W*0v* 

U*2 

u  ^ 

5 

(('onveclive  nncl  non-ccniveciive  il.nta  not  nveragccl) 

4 

.7500 

0.5595 

0.7.S94 

3.2271 

~1.7.7.S6 

3 

.5161 

1.1011 

0.2825 

1  ..7275 

-2.4996 

2 

.2.581 

1.0790 

0.7118 

0..7018 

0.4067 

1 

.0403 

1.0000 

1.0000 

-0.9137 

0.4064 

Covariances  (45  - 

60  km) 

Zi  =  1140m,  w*  = 

1.4445  ms 

iC=  0.3814  ms  ' 

,  0^,^^  =  0.0495 

■K,  q'*’  =  0.0.347  g  kg 

Flight 

Z 

\Cq' 

w'0v' 

u'w' 

v'w' 

Level 

Zi 

w*q* 

W*0v* 

U*2 

5'^ 

1.0807 

0.2786 

0.16.35 

0.5560 

-0.5404 

4 

.81.58 

0.6190 

0.8869 

3.8621 

-0.7.728 

3 

..5614 

0.8494 

0.2710 

1.6478 

-1.7.792 

2 

.2807 

0.8469 

0.6815 

-0.0680 

0.1816 

1 

.0439 

1.0000 

1.0000 

-0.9847 

0.1714 

Covariances  (60  - 

75  km) 

Zi=  1 

340m,  w*  = 

1.4754  ms 

' ,  u*=  0.2493  ms  ^ 

',  0,^  =0.0440 

^K,  q=^=0.0194gkg 

Flight 

Z 

w'q' 

w'0v' 

U'W'' 

v'w' 

Level 

Zi 

w'^q'^ 

W*0v* 

u"2 

s 

.9194 

0.1706 

0.4411 

0.3696 

1.2031 

4 

.6940 

0.5878 

0..5616 

-0..S0I2 

3.6597 

3 

.4776 

1.2249 

0.2863 

0.7666 

-3.7191 

2 

.2.\88 

1.4032 

0.72.74 

I  ..7672 

1.9726 

1 

.0373 

1 .0000 

1.0000 

0.4311 

0.9007 

*  ]  c'vcl  5  is  nhovc  the  convective  IMU  . 

Figttres  3.13  and  3.14  show  the  vertical  profiles  of  the  average  normalized 


crosswind  and  alongwind  stresses  rcspeclivelv  for  the  LESS  and  AMTEX  data.  The 


66 


crosswind  component  generally  increases  with  height  in  the  l.FSS  with  the  exception 
of  the  60-75  km  segment  and  level  5  for  the  45-60km  segment.  Results  for  the 
two  studies  are  in  good  agreement  below  about  0.4  Z/Z|,  indicating  that  crosswind 
shear  is  comparable  at  low  levels  lor  the  two  studies.  Uouever.  the  LESS  crossu  ind 
stress  is  significantly  stronger  from  0.5  Z/Zj  to  approximately  0.9  Z/Z,.  In  order 
to  physically  interpret  the  meaning  of  the  average  covariances,  each  covariance  sum 
was  partitioned  into  its  positive  and  negative  part.  Naturally,  the  one  with  the  larger 
absolute  value  determined  the  sign  of  the  covariance.  1  hen  the  percentage  of  each 
partition  that  had  positive  w'  and  negative  w''  values  was  calculated.  In  this  way, 
the  overall  sign  of  the  other  variable  (u'  or  v')  could  also  be  determined.  From  this 
analysis  it  was  then  possible  to  physically  interpret  the  meaning  of  the  average  covaria¬ 
nces  in  terms  of  upward  or  downward  transport  of  casterly/westerly  (u'<0.  u'>0) 
or  northerly/southerly  (v'<0,  v’>0)  momentum.  Some  of  the  results  of  this  analysis 
and  corresponding  physical  interpretation  are  shown  in  Table  3.3  following  the  verti¬ 
cal  profiles.  In  Table  3.3  the  interpretation  for  easterly  (westerly)  momentum  trans¬ 
fer  refers  to  air  moving  from  east  to  west  (west  to  east).  Similarly,  northerly  (souther¬ 
ly)  momentum  transfer  refers  to  air  moving  from  north  to  south  (south  to  north). 
The  temptation  to  visualize  structure  based  on  this  analysis;  should  be  avoided  since 
results  represent  averages  over  long  flight  segments. 

For  the  normalized  alongwind  component  (v'w'/iF^)  weak  positive  values  are 
found  below  about  0.3  Z/Zi  for  the  LESS  with  strong  negative  values  dominating 
in  the  middle  portion  of  the  PBL  from  0.3  to  0.8  Z/Zj.  Results  for  the  LESS  are 
again  in  good  agreement  at  low  levels  with  those  of  the  AM'I  FX  indicating  that  shears 
were  comparable  (but  slightly  stronger  for  the  AMTFX)  near  the  surface  for  the  two 
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Figure  3.14  Average  normalized  alongvvind  compo?icnl  of 
stress  for  (he  LESS  and  AMTEX  data. 


70 


v'w' 


Figure  3.16  Average  normali/ecl  alongwincl  sliess  for  Ihe 
LESS  and  MASEX  data.' 
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Tahic  3.4  RcsuKs  from  analy  sis  of  sign  of  u',  v'  and  w'  for  covariances. 


Covariances  (45  -  75  km) 
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Sign  of 

Percentage  of 

11 'w'  with 
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1 

+ 
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of  u'/  w' 
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5 

+ 
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4 

+ 

36.41  63.59 

•  /  -( 
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3 

+ 
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momentum 
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Inlet  prelation 

5 

3  0  57  60  !  3 

-  •  k 
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4 
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Table  3.4,  conlinucd 


Percentage  of 
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V  ’  vv ' 

vvitli 
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4 
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ihc  LESS  data  may  he  due  to  the  secondary  maximum  iti  huoyancy  at  this  level  (to 
be  discussed  later).  The  very  strong  alongwind  stress  from  level  2  to  level  3  (Figure 
3.16)  may  be  due  to  the  transition  from  the  suhcloud  to  cloud  layer  at  these  levels. 

Stull  (1988)  points  out  that  buoyancy  is  one  of  the  driving  forces  for  turbu¬ 
lence  in  the  PEL  and  results  from  thermals  of  warm  air  rising  because  they  arc  less 
dense  (positively  buoyant)  than  there  surrounding  environment.  Buoyancy  is  general¬ 
ly  expressed  in  the  following  way: 


gZj  (w'0'v) 


where  g  w'0'v  /  0v  is  called  the  buoyancy  flux.  To  see  why  this  is  equivalent  to 
the  normalized  covariance  of  w  and  ©v,  (i.e.  w'0'v  /  w*  0*v).  substitute  the  expres¬ 
sion  for  w*  (equation  3.1),  the  free  convective  scaling  velocity,  into  (3.8).  After 
simplification,  (3.8)  reduces  to: 


w'0'v  _  w'0'v 

(\7^v)l  W*®' 


(3.9) 


The  profiles  of  normalized  buoyancy  for  the  AMTEX  and  the  LESS  are  in 
very  good  agreement  in  the  lower  half  of  the  PEL.  with  nearly  a  linear  decrease  in 
the  buoyancy  with  height  up  to  0.5  Z/Zj  (Figure  3. 17)  for  both  observational  cases. 
In  the  case  of  the  AMTEX.  the  normalized  buoyanev  continues  to  decrease  from 
this  height  upward,  but  it  does  so  more  slowly  than  in  the  lower  half  of  the  PEL. 
For  the  LESS  results,  however,  a  sharp  increase  in  the  buovanev  is  seen  at  level 
4  for  all  segments  but  espeeiallv  so  for  the  45  -  75  km  and  45  -  60  km  segments. 
The  relatively  large  normalized  buovanev  value  at  level  4  (nearly  90'"r  of  the  surface 


Project  LESS 
'  45  -  75  km 
!  j  45  -  6(i  km 
A  60  -  75  km 


AMTEX 
■  Feb  15  1975 
▼  Feb  16  1975 
•  Feb  18  1975 
A  Feb  22  1975 


Nonnali/.ccl  Buoyancy 

Figure  3.17  Vertical  profile  of  the  average  normalized 

buoyancy  for  the  AMTEX  and  the  LESS  field 
programs. 

value  for  the  45  -  bOkm  segment)  is  most  likely  clue  to  increased  buoyancy  caused 
by  the  release  of  latent  heat  clue  to  condensation  resulting  from  formation  of  precipita¬ 
tion  within  the  clouds. 
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w*©%. 

Figure  3.18  Vertical  profile  of  the  average  nonnalizccl  buoyancy 
for  the  MASEX  and  the  LESS  field  programs 
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Vertical  profiles  of  the  average  normalized  vertical  edclv  hnoyancy  flux  for  both  the 
MASEX  and  LESS  field  programs  are  shown  in  Figure  3. 1  <S.  Again,  good  agreement 
is  generally  found  between  the  two  studies  in  the  lower  portion  of  the  PBL.  The 
data  point  for  the  MASEX  for  AB(W)  at  1.0  Z/Z,  is  dismissed  by  Chou  et  al.  as 
coming  from  a  biased  sample.  According  to  them  this  sample  included  two  cold 
thermals  which  were  associated  with  an  cxircmelv  large  positive  heat  flux,  due  to 
descending  return  flow  of  a  substantial  amount  of  unmixed  air.  By  removing  the 
contribution  of  these  two  thermals,  the  heat  flux  was  reduced  from  90^  to  of 
the  surface  value  (Chou  ct  al..  1986).  With  the  exception  of  this  data  point  the 
normalized  buoyancy  generally  continues  to  decrease  with  height  above  0.5  Z/Zj. 

The  vertical  profiles  of  the  average  normalized  vertical  eddy  moisture  flux 
from  the  MASEX  and  the  LESS  are  shown  in  Figure  3.19.  The  profiles  are  very 
similar  for  the  two  case  studies  with  the  moisture  flux  generally  decreasing  with  height. 
The  especially  strong  fluxes  at  level  2  and,  to  a  lesser  extent,  level  3  for  the  60-75 
km  segment  are  due  to  large  positive  turbulent  q  values  at  these  two  levels  as  evi¬ 
denced  by  the  larger  q  variances  for  this  segment  at  these  levels. 


3.6.2  Comparison  of  LESS  Covariances  with  Model  Results 

A  comparison  between  the  vertical  profile  of  the  average  normalized  buoyan¬ 
cy  for  the  LESS  and  Moeng's  l.ES  model  is  showm  in  Figure  3.20.  Both  the  model 
and  the  LESS  results  are  in  quite  good  agreement  in  the  lower  half  of  the  PBL  hut 
do  not  agree  very  well  in  the  upper  portion.  As  mentiemed  previously,  the  LESS 
profile  shows  a  secondary  maxima  in  the  buoyancy  at  level  4  which  is  due  to  conden 
sation-induced  buovancy.  The  LES  model  profile  does  not  reflect  this  secondary 
maximum,  probably  because  (he  moisture  field  effects  were  excluded  from  the  mod¬ 
el. 
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Figure  3.21  shows  ihe  vcrlical  profile  of  the  non-normnli/cd  huoyanev  for  ihc  Dcar- 
ciorff  model  and  for  the  Project  LESS  60-75km  segment  (note  the  ordinate  is  this 
Figure  is  geometric  height).  The  two  profiles  arc  \erv  similar  with  a  primary  maxi¬ 
mum  near  the  surface  where  healing  is  the  greatest,  followed  hy  a  decrease  with  height 
extending  from  the  subcloud  layer  to  the  cloud  layer.  Both  profiles  indicate  a  second¬ 
ary  maximum  which  is  weaker  than  its  rcspccti\'c  surface  counterpart  (but  still  signifi¬ 
cant)  that  is  most  likely  due  to  the  release  of  the  latent  heal  of  condensation  in  the 
upper  portion  of  the  PBL. 


Buoyancy  x  10-  nrs 


Vertical  profile  of  the  average  non-normali/ed  buoyancy 
for  Dca.dorff's  model  and  the  Project  LHSS  60  75km 
segment. 
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3.7  Turbulence  Kinetic  Energy 

No  analysis  of  a  convective  PEL  would  be  complete  without  a  discussion  of 
the  turbulent  kinetic  energy  (TKE).  The  TKE  provides  a  measure  of  the  intensity 
of  turbulence  and  is  directly  related  to  the  momentum,  heat,  and  moisture  transport 
through  the  boundary  layer  (Stull.  1988). 

The  turbulence  kinetic  energy  per  unit  mass  (TKE/m)  is  defined  in  Ec]uation 
3.10  below: 

TKE  _  Jl_  (u'^  +  v'^  +  w'^) .  (3.10) 

m  2 

For  the  sake  of  brevity,  the  turbulence  kinetic  energy  per  unit  mass  will  be  hereafter 
referred  to  as  simply  TKE  with  (he  “per  unit  mass"  understood.  It  isn't  difficult 
to  sec  the  relationship  between  the  TKE  and  the  variances  of  the  u,  v  and  w-wind 
components  shown  earlier  in  this  chapter;  the  TKE  is  equal  to  one-half  the  sum 
of  these  variances.  The  average  TKE  then  is  just  one-half  the  average  of  the  sum 
of  variances  of  the  three  wind  components  as  shown  in  equation  3.1  1: 


TKE  =  _1  (u'^  +  v'^  +  w'2)  (3  1  1) 

2 

In  a  Type  I  CTBL,  turbulent  kinetic  energy  is  produced  bv  buovanev  and/or  wind 
shear  and  is  transported  by  the  mean  wind  (advection)  and/or  bv  the  turbulent  eddies 
(turbulent  transport).  TKE  may  also  be  redistributed  bv  the  pressure  perturbations 
taking  place  within  the  PBL,  TKE  is  destroyed  or  lost  through  viscous  dissipation 
by  the  smallest  eddies  through  the  energy  cascade  pror  ess  described  earlier. 

The  TKE  can  be  easily  broken  down  inti>  its  components  and  other  terms 
to  further  illustrate  the  makeup  of  IKE  as  described  in  the  equations  below: 
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Average  horizontal  eomponent  of  TKE:  TKE|,  =  1  (ii'^  +  v'^)  (3.12) 


Average  vertical  component  of  TKE:  TKEv  =  1  (3.13) 


Averase  vertical  eddv  w'TKE  =  w'  (u'^  +  \  +  \\'^)  (3.14) 

2 

w'TKEh  =  v^ )  (3.15) 


transport  of  TKE. 

Average  vertical  eddy  transport 
of  the  horizontal  component  of 
TKE: 


Average  vertical  eddy  transport  w'TKE,  =  vr  ixT-)  (3.16) 

of  the  vertical  eomponent  of  “ 

TKE: 


Ratio  of  the  average  vertical 
eddy  transport  of  the  horizontal 
component  of  TKE  to  the  aver¬ 
age  vertical  eddy  transport  of 
the  vertical  component  of  TKE: 


w'TKEi,  =  \v'“  (u'^+T^  (3.17) 

w'TKEv  w'  (u'^ 


Table  3. .5  Project  LESS  boiiiKlar\'  layer  normalized  average  TKE  statistics 


TKE  Statistics  (45  -  75  km) 


Zi  = 
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'  =  l.46<S4  ms  ' 

z 

TKE 

TKEh 

1KE„  wTKE 

w'TKEh 

w  TKE^. 

w'TKEh 

Zi 

VV*2 

w*-’ 

w*'" 

w^’ 

w'TKE, 

0.993.S 

(C’nnvcclivc  ;iikI  non 

-convective 

n<'t  jivei 

r.tgeil) 

0.7.S0O 

0,7HR7 

0.,S0,77 

n.zsso 

0.0720 

0,0008 

0.07 1 1 
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Tabic  3.5,  conlinuccl 


TKE  Stalislics  (45  -  bO  km) 
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3.7.1  Comparisons  among  Obscrvalional  Case  Suiciics. 

Figure  3.22  shows  the  vertical  profile  of  the  normalized  average  turbulence 
kinetic  energy  for  the  Protect  l.F.SS.  Maxima  arc  feumd  for  ail  three  segments  m 
level  2  and  at  level  4  for  the  45  -  75  km  and  45  -  60  km  segments.  Intercstingiv, 
the  maximum  at  level  2  does  not  occur  at  the  level  at  which  buovanev  is  a  maximum 
(level  1,  Figure  3.20).  The  vertical  profiles  of  the  average  normalized  horizvmtai 
(Figtirc  3.23)  and  vertical  (Figure  3.24)  components  (4  I  KIv  indicate  that  the  maxi¬ 
mum  at  level  2  is  due  primarilv  to  a  corresponding  maximum  in  the  horizontal 
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component  of  TKE  at  this  same  level.  The  vertical  profile  of  the  vertical  component 
of  TKE  on  the  other  hand  show's  very  little  variability  in  the  w/  variance  except  at 
level  5.  By  examining  the  u  and  v  variances  which  make  up  the  horizontal  component 
of  TKE  it  becomes  evident  that  the  v-component  of  variance  is  largely  responsible 
for  the  maximum  in  TKE  at  level  2  since  the  v  variance  is  between  5^  (60-75  km 
segment)  and  40%  (45-75  km  segment)  larger  than  the  u-component  of  variance. 
The  second  maximum  at  level  4  for  the  LESS  data  again  appears  to  be  due  to  the 
horizontal  rather  than  the  vertical  component  of  TKE  since  the  horizontal  component 
is  63%  (45  -  60km  segment)  to  76%  (45  -  75km  segment)  larger  than  the  vertical 
component.  However,  unlike  the  maximum  at  level  2,  the  maximum  at  level  4  owes 
its  origin  primarily  to  the  u-component  variance  rather  than  the  v-component.  At 
level  4  the  u  variance  is  between  13%  to  48%  larger  than  the  v  variance.  The  u 
variance  at  level  4  is  more  than  two  times  as  strong  as  the  w— component  for  the 
45  -  60km  segment  but  the  two  are  comparable  for  the  45  -  75km  segment. 

A  comparison  between  the  LESS  and  the  MASEX  data  for  the  normalized 
average  horizontal  and  vertical  components  of  TKE  shows  the  values  of  the  MASEX 
to  be  substantially  larger  (about  200%  larger)  in  the  lower  half  of  the  convective 
PBL  but  of  comparable  value  in  the  upper  one-third.  This  would  indicate  that  for 
a  given  surface  buoyancy  flux  (actually  level  1  buoyancy  flux  for  the  LESS)  as  ac¬ 
counted  for  by  the  normalization  by  w*,  the  MASEX  PBL  is  dynamically  more  active 
than  the  LESS. 

Figure  3.25  shows  the  vertical  profile  of  the  normalized  average  vertical  flux 
of  TKE  for  the  LESS  and  MASEX  cases.  With  a  few  exceptions,  the  vertical  flux 
decreases  with  height  for  the  LESS  case  with  the  maximum  found  at  level  2.  Below 
this  maximum,  there  is  more  upward  flux  leaving  the  top  of  any  given  layer  than 


Figure  3.25  Vertical  profile  of 
flux  of  TKE  for  t 
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enters  from  below,  resulting  in  a  net  divergence  or  loss  of  TKE.  Above  this  maxi¬ 
mum,  there  is  a  net  convergence  or  gain  of  TKE.  As  a  result,  some  of  the  the  TKE 
produced  at  the  lowest  levels  where  buoyancy  is  a  maximum  is  transported  to  higher 
levels.  A  second  maximum  is  found  at  level  4  for  the  LESS  profile  for  the  45-60km 
segment.  This  peak  appears  to  be  due  equally  to  the  vertical  flux  of  horizontal  TKE 
for  which  this  .same  maximum  is  found  (Figure  3.26)  and  the  vertical  flux  of  vertical 
TKE  (Figure  3.27).  Results  for  the  LESS  case  study  compare  favorably  with  those 
of  the  MASEX  study  with  both  showing  that  the  vertical  transport  of  the  vertical  com¬ 
ponent  of  TKE  is  two  to  four  times  stronger  than  the  vertical  transport  of  the  horizontal 
component. 

The  vertical  profiles  of  the  normalized  average  vertical  flux  of  the  horizontal 
and  vertical  components  of  TKE  are  shown  in  Figure  3.28  and  3.29  respectively 
for  the  LESS  and  the  AMTEX  results  of  Lenschow  et  al.  (1980).  The  vertical  flux 
profiles  for  the  horizontal  component  of  TKE  are  very  similar  for  the  two  case  studies 
with  the  exception  of  the  aforementioned  strong  maximum  in  the  LESS  at  level  4 
for  the  45  -  60km  segment.  For  the  vertical  flux  of  the  vertical  component  of  TKE, 
however,  the  LESS  values  are,  with  very  few  exceptions,  larger  than  those  for  the 
Lenschow  results  for  the  AMTEX.  This  indicates  that  for  a  given  buoyancy  produc¬ 
tion  the  turbulent  eddies  for  the  LESS  case  arc  more  efficient  and  capable  at  trans¬ 
porting  the  vertical  component  of  TKE  than  arc  those  fo.  the  AMTEX  case.  With 
regards  to  the  vertical  flux  of  horizontal  TKE,  the  turbulent  eddies  appear  to  be  ap¬ 
proximately  equal  for  both  cases  in  terms  of  their  ability  to  transport  the  TKE. 

Figures  3.30  and  3.31  show  the  ratio  of  the  vertical  flux  of  velocity  variance 
between  the  horizontal  and  the  vertical  components  of  TKE  for  the  LESS  and  AM¬ 
TEX  (Figure  3.30)  and  the  LESS  and  MASEX  (Figure  3.31).  For  both  the  AMTEX 
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and  MASEX  cases  the  ratio  is  nearly  unity  at  the  lowest  level  indicating  that  each 
component  contributes  equally  to  turbulent  transport  at  this  level.  For  the  LESS, 
however,  the  ratio  is  less  than  unity  indicating  that  w'-  contributes  more  to  turbulent 
transport  at  this  level  than  for  either  the  AMTEX  or  M,A,SEX  case.  This  is  a  result 
of  a  stronger  vertical  transport  of  vertical  variance  in  the  LESS  over  that  of  the  AM 
TEX  and  MASEX.  The  ratio  drops  to  less  than  unity  in  the  mid-portion  of  the 
PEL  for  both  the  AMTEX  and  the  MASEX  results  indicating  that  ^  contributes 
more  to  turbulent  transport  in  this  portion  of  the  PEL  for  both  cases.  For  the  LESS 
case  the  ratio  remains  nearly  constant  with  height  until  level  5  is  reached  where  the 
ratio  for  both  the  45  -  60km  and  60  -  75km  segments  becomes  negative. 

3.7.2  Comparisons  between  the  LESS  and  Model  results 

Figure  3.32  gives  the  vertical  profile  of  the  normalized  average  horizontal 
component  of  TKE  for  the  Project  LESS  study  and  the  3~D  model  of  Deardorff 
(1980).  The  horizontal  component  is  generally  stronger  for  the  LESS  study  than 
for  Deardorff’s  simulation  of  day  33  of  the  Wangara  experiment  indicating  the  for¬ 
mer  has  stronger  turbulence  intensity  in  the  horizontal  component. 
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3 . 8  Skewness 

In  statistics,  skewness  is  often  used  as  a  measure  of  the  extent  to  which  items 
are  symmetrically  distributed.  One  way  to  study  the  symmetry  (or  asymmetry)  of 
a  frequency  distribution  is  by  examining  the  mode,  median  and  mean  where  the  mode 
is  the  location  on  the  scale  which  has  the  greatest  concentration  of  items  and  the 
median  has  roughly  half  of  the  items  lying  below  and  above  it.  (Neter  et  al.,  1988) 


(a)  symetrical  (b)  skewed  left  (negatively) 


Median 


Median 


Figure  3.33  Examples  of  .symmetrical  and  skewed  imimodal  frequency 
distributions.  From  Neter  et  al.  .  1988. 
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The  three  frequency  polygons  in  Figure  3.33  show  the  rclalivc  positions  of  the  mean, 
median  and  mode  for  a  unimodal  frequency  distribution  with  varying  degrees  and 
direction  of  symmetry.  Figure  3.33a  shows  a  symmetrical  distribution  for  which 
the  mean,  median  and  mode  are  identical.  Figures  3.33b  and  3.33c  show  frequency 
distributions  skewed  to  the  left  (negatively)  and  skewed  to  the  right  (positively)  respec¬ 
tively.  The  mean  in  each  ease  is  typically  furthest  out  toward  the  tail  of  the  distribu¬ 
tion  because  of  its  sensitivity  to  large  and  small  values  and  the  median  is  between 
the  mean  and  the  mode. 

Mathematically,  the  skewness  can  be  expressed  for  a  set  of  items  Xj  for  i  = 
1  to  n  as  shown  below: 


Sx  =  m3 
s 


(3.12) 


where  m3 


(Xj  -  X)-^ 
n-1 


and 


s  = 


N 


n-1 


i  =  I 


In  other  words,  the  standardized  skewness  measure,  Sx,  is  the  average  third  moment 
about  the  mean  divided  by  the  standard  deviation  cubed.  Because  of  the  cube  in 
the  numerator  of  Sx,  large  deviations  Xj  -  X  will  dominate  the  sum  in  the  numerator 
of  Sx.  If  the  large  deviations  arc  mostly  positive,  the  skewness  will  be  positive,  if 
they  arc  largely  negative  the  skcwme.ss  will  be  negative  (the  denominator  is  always 
positive).  Because  large  deviations  appear  in  the  tail  portion  of  the  frequency  distribu¬ 
tion,  one  can  see  from  Figure  3.33  that  skewness  will  he  positive  or  negative  accord¬ 
ing  to  whether  the  direction  of  skewness  is  positive  (right)  or  negative  (left)  (Neter 
et  al.,  1988). 

Skewness  is  a  useful  measure  in  the  PBL  because  it  gives  an  indication  of 
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the  structure  of  the  motion.  For  instance,  if  w  =  0  and  if  Sw  >  0,  updrafts  are 
narrower  and  stronger  than  surrounding  downdrafts  (Moeng  and  Rotunno,  1990). 
In  the  sections  that  follow,  vertical  profiles  of  the  skewness  of  all  five  meteorological 
parameters  (u,  v,  w,  0v  and  q)  will  be  presented  with  comparisons  made  between 
the  LESS  results  and  other  ease  siudics  and  between  the  LESS  and  model  results. 

з. 8.1  U,  V,  0v  and  Q  Skewness 

In  this  section,  vertical  profiles  of  skewness  will  be  shown  and  discussed  for 

и,  V,  0v  and  q  for  the  Project  LESS.  Results  for  the  veriical  velocity,  w,  will  be 
presented  in  this  section,  but  a  discu.ssion  of  the  vertical  profiles  will  be  held  off 
until  the  following  section. 

Table  3.6  summarizes  the  skewness  results  for  the  Project  LESS  data  sets. 
As  has  been  the  case  for  all  the  results  presented  thus  far,  skewness  calculations 
have  been  done  using  data  detrended  using  a  .second-degree  polynomial  fit.  As  a 
result,  the  mean  for  each  variable  is  zero  per  the  earlier  discu.ssion. 

Table  3.6  Project  LESS  boundary  layer  skewnc,ss  values. 

Skewness  (45  -  75  km) 


Zi  =  124()m 
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0.6703 

0.8875 
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Table  3.6,  continued 


Skewness  (45  -  60  km) 

Zi  =  1140m 
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Skewness  (60  - 

75  km) 
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*  I  .evel  5  is  above  the  convective  PFiL. 


The  vertical  profile  of  the  skewness  of  the  u-wincl  component  is  shown  in 
Figure  3.34  below.  Skewness  is  weakly  negative  for  level  1,  indicating  narrow  re¬ 
gions  of  relatively  strong  westerly  momentum  surrotmded  by  broad  regions  of  general¬ 
ly  weak  easterly  momentum.  At  level  2,  just  the  opposite  occurs;  skewness  is  negative 
for  all  three  segments,  indicative  of  relatively  strong  narrow  regions  of  easterly  mo¬ 
mentum  surrounded  by  regions  of  broad,  relatively  weak  westerly  momentum.  With 
the  exception  of  the  60  -  75km  segment  at  levels  2  and  3,  the  u-component  of 
skewness  is  generally  positive  throughout  the  rest  of  the  PBL  and  generally  increases 
in  strength  with  increasing  height. 
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Figure  3.34  Vertical  profile  of  the  u-wind  component  skewness 
for  the  Project  LESS. 
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The  vertical  profile  of  the  v-component  of  skewness  (Figure  3.35)  shows  values  very 
nearly  zero  in  the  lower  half  of  the  PBL  and  positive  values  at  level  4  indicating 
narrow  regions  of  relatively  strong  northerly  momentum  coupled  with  generally  broad 
regions  of  weak  southerly  momentum.  The  skewness  at  level  5  for  both  the  45  - 
60km  and  60  -  75km  segments  is  quite  strongly  negative,  corresponding  to  narrow 
(broad)  regions  of  strong  (weak)  southerly  (northerly)  momentum. 

The  skewness  of  the  virtual  potential  temperature  (Figur?  3.36)  is  strongly 
positive  near  the  surface  and  then  decreases  nearly  linearly  with  height  throughout 
the  remainder  of  the  PBL.  The  strong  positive  skewness  at  the  lowest  level  is  indicative 
of  relatively  narrow  (broad)  regions  of  relatively  warmer  (cooler)  air,  reflective  of 
the  strong  buoyancy  at  this  level. 

The  skewness  of  the  specific  humidity  (Figure  3.37)  is  also  strongly  positive 
at  lower  levels,  corresponding  to  the  strong  moisture  flux  near  the  Lake  Michigan 
surface.  At  the  mid-portion  of  the  PBL,  however,  the  skewness  is  strongly  negative 
indicating  that  there  are  relatively  narrow  regions  of  drier  air  surrounded  by  relative¬ 
ly  broad  regions  of  moi.st  air  within  the  top  of  the  cloud  layer. 

3.8.2  Vertical  Velocity  Skewness 

As  pointed  out  by  Moeng  and  Rotunno  (1990),  a  number  of  puzzling  features 
concerning  vertical  velocity  skewness  have  emerged  from  the  comparison  of  LES 
model  results  with  those  from  observational  case  studies  involving  convective  PBLs. 
Figure  3.38  gives  the  vertical  profile  of  vertical  velocity  skewness  for  the  LESS,  AM- 
TEX  and  LES  results  of  Moeng.  For  the  LESS  case,  skewness  is  a  maximum  at 
level  1,  decreases  slightly  with  height  within  the  subclouci  layer  and  then  is  nearly 
constant  with  height  in  the  cloud  layer. 
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Figure  3.38  Vertical  profiles  of  the  vertical  velocity  skewness  for  the  LESS 
and  the  AMTEX  (solid  circles)  and  Moeng's  LES  model 
(heavy  dark  line)  for  surface  healing  only. 
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The  solid  circles  for  the  AMTEX  results  are  taken  from  Moeng  and  Rotunno  (1990) 
and  agree  well  with  the  LES  profile.  The  results  from  Moeng’s  LES  model  (denoted 
by  the  heavy  dark  line)  are  for  the  case  of  surface  heating  only.  Skewness  values 
for  the  LES  model  are  smaller  than  those  of  the  AMTEX  or  LESS  study  near  the 
surface  probably  because  subgrid-scale  features  dominate  the  LES  results  at  these 
levels.  The  skewness  for  the  model  continues  to  increase  until  it  reaches  a  maximum 
at  approximately  0.9  Z/Zj.  It  is  in  the  upper  portion  of  the  PEL  that  model  and 
observational  case  results  generally  do  not  agree.  For  example,  as  shown  in  Figure 
3.38,  the  skewness  remains  nearly  constant  for  the  AMTEX  and  the  LESS  results 
in  the  upper  half  of  the  PEL  but  for  Moeng's  model  continues  to  increase.  This 
increase  in  Sw(7.)  exits  in  all  other  large  eddy  simulations  of  the  convective  PEL 
(Moeng  and  Rotunno,  1990).  This  discrepancy  in  the  upper  portion  of  the  PEL 
prompted  Moeng  and  Rotunno  to  further  investigate  the  cause  of  the  increase  in  posi¬ 
tive  skewness  with  height  in  the  upper  portion  of  the  PEL  for  a  layer  of  thermal 
convection.  They  conducted  numerical  simulations  for  both  no-slip  and  frcc-slip 
top  conditions  for  Ra=380,000  and  Pr=l.  Surprisingly,  they  found  that  the  no-slip 
condition  results  most  closely  replicated  the  LES  simulation  results  for  surface  heating 
only.  The  free-slip  condition  for  the  surface  heating  case  resulted  in  a  skewness 
profile  that  is  positive  and  relatively  constant  above  a  certain  height.  This  surprising 
result  prompted  Moeng  and  Rotunno  to  further  investigate  whether  the  LES  model 
contains  some  extra  subgrid-scale  stress  near  the  upper  part  of  the  convective  zone 
that  would  make  its  skewness  behave  more  like  that  of  the  no-slip  rather  than  the 
free-slip  condition.  Their  examination  of  the  data  led  to  an  answer  in  the  negative. 
They  went  on  to  show  why  they  believed  the  temperature  inversion  typically  overlaying 
a  convective  PEL  alters  in  some  way  the  convective  eddies  impinging  on  this  layer 
in  a  way  that  resembles  the  no-sIip  experiment.  To  test  this  hypothesis  Moeng  and 
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Rotunno  carried  out  a  third  experiment  with  surface  heating  only  with  an  imposed 
inversion  in  the  upper  portion  of  the  numerical  domain.  The  resulting  vertical  skew¬ 
ness  profile  (Figure  1 1 ,  after  Moeng  and  Rotunno,  1 990)  is  very  similar  to  the  large- 
eddy  simulations  of  a  surface-heating-driven  PBl-  and  a  surfacc-hcating-drivcn 
PBL  with  a  no-slip  top.  They  explain  the  increase  in  skewness  at  the  top  of  the 
PBL  by  considering  the  effect  of  the  inversion  on  the  turbulent  kinetic  energy.  They 
hypothesize  that  the  updrafts  lose  their  kinetic  energy  at  the  inversion  resulting  in 
only  broad  weak  downdrafts  and,  hence,  positive  skewness.  Mason  (1989)  hypothe 
sized  that  observed  skewnesses  at  the  top  of  the  mixed  layer  are  lower  than  the  simu 
lated  ones  due  to  the  presence  of  “large  scale  features  and  inhomogeneitics"  in  the 
real  atmosphere.  Mason  went  on  to  say  that  skewness  in  the  boundary  layer  can 
be  lowered  below  that  replicated  in  simulations  by  the  presence  of  mcsoscalc  features 
if  (a)  these  features  are  sinusoidal  and  therefore  contribute  to  total  vertical  velocity 
variance  but  have  near-zero  skewness;  and  (b)  if  these  features  cannot  be  replicated 
by  the  simulation  in  question.  These  features  may  be  related  to  gravity  waves  interact¬ 
ing  with  the  boundary  layer  as  described  by  Clark  el  al.  (1984).  There  is  some 
evidence  of  an  increase  in  skewness  at  level  4  for  the  LESS  45-60km  segment  but 
this  is  most  likely  due  to  the  secondary  maximum  in  buovancy  at  this  flight  level. 


4.  DETRENDED  VS.  RAW  DATA 


As  mentioned  earlier  in  Chapter  2,  Wilc/.ak  and  Businger  (1983)  have  shown 
that  when  computing  lurhulencc  statistics  in  the  convective  PBL  there  is  little  differ¬ 
ence  between  linear  and  quadratic  detrending  when  compared  to  the  variability  when 
no  detrending  is  used.  I  his  relatively  short  chapter  will  examine  vertical  profiles  of 
a  variety  of  turbulent  statistics  for  which  both  second-degree  polynomial  detrended 
and  non-detrended  data  have  been  used.  In  the  latter  case,  equation  (3.5)  still  holds 
where  the  average  subtracted  from  the  instantaneous  value  of  the  variable  is  simply 
the  Reynolds  average  rather  than  the  predicted  value  for  the  given  meteorological 
variable. 

4.1  TKE 

Figure  4. 1  shows  the  vertical  profile  of  the  normali/cd  average  TKE  for  both 
the  detrended  (open  symbols)  and  non-detrended  or  raw  (closed  symbols)  data  of 
the  Project  LESS.  The  two  sets  of  curves  arc  both  qualitativclv  and  quantitatively 
in  good  agreement,  with  the  TKE  values  for  the  non-detrended  data  consistently 
slightly  larger  than  those  for  the  detrended.  This  certainly  is  not  surprising  since  one 
would  expect  the  second  moments  calculated  about  the  mean  (non-detrended)  to 
be  larger  than  those  calculated  about  the  detrend  curve.  Generally  maxima  and 
minima  arc  found  at  the  same  levels  in  both  the  detrended  and  raw  data  sets.  The 
fact  that  the  two  sets  of  curves  are  quantitatively  in  good  agreement  suggests  that  very 
little  synoptic  or  largc-mesoscale  signal  is  present  in  the  u.v  and  u-wind  components. 


4.2  Normalized  Buoyancy 

The  vertical  profiles  of  the  normalized  buoyancy  for  the  detrended  and  raw 
data  (Figure  4.2)  present  a  picture  quite  different  from  that  of  the  TKE.  The  iuioyan- 
cy  profiles  for  the  detrended  data  for  all  three  segments  arc  remarkably  similar  with 
the  aforementioned  maximum  found  at  levels  1  and  4.  The  striking  similarity  among 
the  curves  for  the  detrended  data  indicates  that  after  synoptic  and  largc-mesoscale 
signals  are  removed  from  the  data,  the  45-75Km  “hand-picked”  data  segment  repre¬ 
sents  a  nearly  homogeneous  portion  of  the  convective  PBL.  The  raw  data,  on  the 
other  hand  do  not  indicate  the  same  convergence  in  the  buoyancy  statistics.  In  fact, 
the  two- 15  km  segments  arc  nearly  mirror  images  of  one  another  with  maximum 
in  the  normalized  buoyancy  for  the  45-60km  segment  (at  level  2)  occurring  at  the 
same  level  as  the  minimum  for  the  60-75km  segment  for  flight  levels  below  0.8 
Z/Zj.  Likewise,  the  maximum  in  normalized  buoyancy  for  the  60-75km  segment 
occurs  at  the  same  level  as  the  minimum  for  the  45-fiOkm  segment  buoyancy  for 
levels  below  0.8  Z/Zj.  The  lack  of  convergence  in  the  normalized  buoyancy  statistic 
in  this  case  suggests  that  synoptic  and  largc-mesoscale  features  may  be  affecting  the 
b  oyancy.  Considering  the  lack  of  trend  in  the  vertical  velocity,  it  would  seem  that 
the  large  scale  features  arc  more  likely  found  in  the  virtual  potential  temperature. 
The  comparison  of  the  normalized  average  variance  of  virtual  potential  temperature 
for  both  the  raw  and  detrended  data  (Figure  4.3)  seems  to  support  this  hypothesis. 
The  variances  for  the  raw  data  for  virtual  potential  temperature  arc  as  much  as  300'"r 
larger  than  tho.sc  for  the  detrended  data. 

4.3  Vertical  Velocity  Skewness 

Figure  4.4  shows  the  vertical  profile  of  the  vertical  velocity  skewness  for  the 
Project  LESS  detrended  and  raw  data.  As  was  the  ca'^c  for  TKE.  the  two  sets  of 
curves  arc  in  very  good  agreement  both  qualitatively  and  quantitatively. 


z, 


Figure  4.2  Vertical  profile  of  the  normali/ed  buoyancy  for  both 
(Jctreiided  and  non-detrended  data  for  Project  LESS. 
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Figure  4.4  Vcrlical  profile  of  the  vertical  velocity  skewness  for  the  Pro¬ 
ject  LESS  for  both  the  detrended  and  raw  data. 


As  pointed  by  out  by  Agee  and  Hart  (1990),  the  vertical  velocity  skewness  can  be 
easily  expanded  into  three  terms  as  shown  below: 


=_^Z_  +  +  2w~ 

(^)3/2  (4.1) 

or  S  =  S,  +  S2  +  S3  (4-2) 

Because  of  continuity  constraints  numerical  models  have  the  condition  that  w=0, 
which  effectively  reduces  (4.2)  to  S  =  S|.  This  is  also  the  case  for  the  detrended 
vertical  velocity  data  for  which  w  =  0  per  the  earlier  discussion.  In  the  real  atmospher¬ 
ic  convective  boundary  layer  this  constraint  of  course  docs  not  exist.  It  seems  some¬ 
what  surprising  then  that  the  skewness  profiles  for  the  raw  and  detrended  vertical 
velocity  data  are  so  similar.  Table  4.1  shown  below  gives  the  relative  contributions 
of  Si,  S2  and  Sato  the  total  skewness  for  all  three  flight  segments  of  the  LESS: 


Table  4.1  Relative  contributions  of  the  S|,  S2  and  S3  terms  to  (he  total  skewness 
of  vertical  velocity. 


45  -  75  km  Segment 


Z/Zi 

W 

Si 

S2 

S3 

S 

0.7500 

0.1182 

0.7395 

-0.3141 

0.0022 

0.4277 

0.5161 

0.3192 

1.2038 

-0.9037 

0.0435 

0.3437 

0.2581 

0.2035 

1.2385 

-0.5125 

0.0092 

0.7351 

0.0403 

-0.0220 

0.7962 

0.0608 

0.0000 

0.8569 

45  -  60  km  Segment 


Z/Zi 

W 

s, 

S2 

S3 

S 

1.0807 

-0.0307 

2.9555 

0.1379 

-0.0002 

3.0933 

0.8158 

0. 1900 

0.8834 

-0.4807 

0.0077 

0.4103 

0.5614 

0.1946 

0.9320 

0.5538 

0.01  14 

0.3896 

0.2807 

0.1898 

1.2483 

-0.4808 

0.0077 

0.7751 

0.0439 

-0.0419 

0.7638 

0. 1134 

-0.0001 

0.8771 

Tabic  4.1,  continued 

60  -  75  km  Segment 


Z/Zi 

W 

s, 

S2 

S3 

S 

0.9194 

-0.1836 

-0.6658 

1.0906 

-0.0707 

0.3541 

0.6940 

0.0461 

0.5180 

-0. 1310 

0.0002 

0.3872 

0.4776 

0.4353 

1.4800 

-1.3050 

0.1098 

0.2849 

0.2388 

0.2172 

1.2262 

-0.5429 

0.0109 

0.6942 

0.0373 

-0.0022 

0.8327 

0.0061 

0.0000 

0.8388 

Generally,  the  Si,  which  is  proportional  to  the  tuibulcnt  vertical  transport  of 
vertical  velocity  variance,  makes  the  largest  contribution  to  the  total  skewness  although 
the  contribution  of  the  Sa  is  occasionally  comparable  to  and  may  even  exceed  the 
Si  term  contribution.  Because  the  S3  involves  the  cube  of  the  average  vertical  velocity 
(which  is  always  less  than  unity)  its  contribution  to  total  skewness  is  the  lea.st  signifi¬ 
cant.  As  pointed  out  by  Agee  and  Hart  (1990),  however,  the  S2  and  S3  terms  may 
make  important  contributions  to  skewness  when  w  values  exceed  0.1  ms“'  for  obser¬ 
vational  periods  <10  minutes  (the  longest  observational  period  for  the  LESS  is  less 
than  8  minutes).  This  is  the  case  because  the  vertical  velocity  data  meet  specific 
accuracy  requirements  (see  NCAR  RAF  Bulletin  23).  It  would  appear  from  this 
analysis  that  the  vertical  velocity  skewness  for  raw  and  detrended  data  arc  very  similar 
because  of  the  lack  of  significant  trend  in  this  component  of  the  wind.  It  appears 
that  Si  for  the  detrended  data  effectively  equals  S|  +  S2  +  S3  for  the  raw  data  and, 
hence,  the  total  skewnesses  arc  nearly  the  same  when  this  condition  is  met. 

This  section  is  brought  to  a  close  with  a  closer  analysis  of  the  sign  and  magni¬ 
tude  of  the  vertical  velocity  for  the  Project  LESS  data.  Table  4.2  gives  an  account 
of  the  overall  average  vertical  velocity,  w  (always  zero  for  detrended  data),  the  num¬ 
ber  of  positive  and  negative  vertical  velocities,  the  average  of  positive  and  negative 
Wj  values  and  the  summation  of  the  length  segments  for  the  positive  and  negative 
values. 


Table  4.2  Results  from  analysis  of  the  sign  and  inagniliidc  of  the  vertical  velocity 
for  both  detrended  and  raw  data  of  the  Project  LESS  for  the  45  -  75  km  segment. 


Detrended  Data 


45 

-  75  km  Segment 

Total 

#of 

#  of 

Average 

Average 

I  ,ength  of 

Length 

Z/Zi 

#  obs 

pos.  w 

neg.  w 

of  pos.  w 

of  ncg.vv 

pos.  w(m) 

neg.  w( 

0.7500 

8441 

4053 

4388 

0.9029 

-0.8344 

14405 

15595 

0.5161 

8063 

3735 

4328 

0.9416 

-0.8123 

13897 

16103 

0.2581 

8256 

3586 

4688 

1.0702 

-0.8147 

12965 

17035 

0.0403 

8683 

3684 

4999 

0.9811 

-0.7230 

12728 

17272 

Raw  Data 

45 

-  75  km  Segment 

Total 

#  of 

#  of 

Average 

Average 

Length  of 

Length 

Z/Zi 

#  obs 

pos.  w 

neg.  w 

of  pos.  w 

of  neg.w 

pos.  w(m) 

neg.  w| 

0.7500 

8441 

4461 

3980 

0.9555 

-0.8199 

15855 

14145 

0.5161 

8063 

4710 

3353 

1.0414 

-0.7093 

17525 

12475 

0.2581 

8256 

4180 

4076 

1.1302 

-0.7471 

15189 

14811 

0.0403 

8683 

3612 

5071 

0.9783 

-0.7346 

12480 

17520 

For  the  detrended  data,  the  analysis  of  the  sign  and  magnitude  of  w  is  consis¬ 


tent  with  the  values  of  the  vertical  velocity  skewness.  For  example,  the  skewness 
at  level  4  for  the  45  -  75  km  .segment  was  found  to  be  0.3305  (Table  3.4).  For 
the  detrended  data,  there  are  fewer  positive  than  negative  vertical  velocities  at  this 
level  and  the  average  (-t-)  vertical  velocity  is  stronger  (-t-0.9029)  than  the  average 
(-)  vertical  velocity  (-0.8344).  The  summation  of  the  length  .segments  for  positive 
and  negative  vertical  velocities  confirms  that  this  level  is  marked  by  relatively  narrow 
regions  of  strong  updraft,  surrounded  by  relatively  broad  regions  of  downdraft.  The 
results  for  the  raw  data  do  not  lead  to  the  same  conclusion.  For  the  same  level  4 
in  the  raw  data,  (+)  velocities  are  more  numerous  but  still  stronger  than  (-)  velocities. 
Hence,  one  would  conclude  that  at  this  flight  level,  broad  regions  of  relatively  stronger 
upward  vertical  velocities  were  dominating  over  narrow  regions  of  relatively  weak 
downdrafts.  This  analysis  points  out  the  importance  of  not  drawing  conclusions  about 


the  structure  of  updrafts  and  downdrafts  from  vertical  velocity  skewness  when  vv  does 
not  equal  zero. 
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5.  SUMMARY  AND  CONCLUSIONS 

Mcleorlogical  data  sampled  at  a  frequency  of  20Hz  and  collected  by  NCAR  research 
aircraft  have  been  analyzed  for  the  CAO  event  of  10  Jan  84  over  Lake  Michigan  during 
the  Project  LESS.  This  particular  CAO  event  was  associated  with  a  surface  anticyclone-cy¬ 
clone  pair  between  which  southward  transport  of  continental  polar  air  occurred  resulting 
in  the  formation  of  a  Type  I  CTBL  over  the  relatively  warm  waters  of  Lake  Michigan. 
A  brief  comparison  of  850mb  temperatures  from  the  PROAM  and  PROAMU  for  the  10 
Jan  84  case  and  the  results  from  the  analysis  done  by  Konrad  and  Colucci  (1989)  indicates 
this  event  represents  a  strong  if  not  extreme  CAO  event. 

The  NCAR  King  and  Queen  Air  aircraft  collected  data  at  five  vertically  stacked  flight 
levels  oriented  perpendicular  to  the  mean  wind  (which  was  northerly)  at  a  latitude  of  43°N 
in  a  west-east  and  cast-west  fashion  (i.c.  perpendicular  to  the  2-d  cloud  streets  which  had 
formed).  Because  it  was  desirable  to  obtain  turbulence  statistics  over  as  homogeneous  a 
region  as  possible,  a  “hand-picked”  data  set  that  attempted  to  eliminate  the  effects  of  the 
interfacial  and  inversion  layer  and  the  Michigan  and  Wisconsin  land  breezes  was  chosen. 
The  resulting  data  set  included  the  u,v,w,0v  and  q  20Hz  data  collected  45-75km  from  the 
Wisconsin  shore.  This  turned  out  to  be  a  fortuitous  choice  since  the  data  collection  period 
over  this  flight  segment  was  less  than  10  minutes  for  each  flight  level,  allowing  for  the  use 
of  mean  vertical  velocities  in  the  Sw  calculations  for  raw  data. 

Data  used  in  the  turbulence  statistics  have  been  detrended  using  a  second-degree 
polynomial  fit.  Comparisons  were  made  between  the  LESS  results  and  those  for  two  other 
observational  case  studies;  the  AMTEX  and  MASEX.  Turbulence  statistics  were  also  com¬ 
pared  to  those  of  the  models  of  Deardorff  (1980)  and  Moeng(1984).  Comparisons  among 


eibscrvational  case  studies  and  model  results  are  facilitated  by  the  use  cif  appropriate  scaling 
factors  such  as  Z,  and  w'^,  the  depth  of  the  mixed  layer  and  the  free  convective  scaling 
velocity,  respectively. 

Vertical  profiles  of  the  means  of  u  and  v  show  very  little  shear  in  the  mean  horizontal 
wind  except  near  the  base  of  the  inversion  at  the  top  of  the  mixed  layer.  The  vertical  profile 
of  the  mean  vertical  velocity  shows  a  positive  maximum  in  the  subcloud  layer  and  weak 
subsidence  near  the  inversion  base  for  the  45-60km  and  60-75km  segments.  The  vertical 
profile  of  the  mean  virtual  potential  temperature  and  specific  humidity  reveals  the  drying 
and  warming  effect  of  the  entrainment  of  air  from  the  inversion  layer  overlying  the  convec¬ 
tive  PBL. 

The  vertical  profiles  of  the  normalized  average  u,v  and  w  variances  indicates  that 
the  values  are,  generally  speaking,  larger  for  the  LESS  results  than  for  the  AMTEX,  especial¬ 
ly  for  the  w  variance.  This  would  indicate  that  for  a  given  surface  buoyancy  flux,  as  ac¬ 
counted  for  in  the  normalization  by  w*.  the  LESS  CTBL  is  more  turbulent  that  the  AMTEX. 
The  vertical  profiles  of  the  normalized  average  variances  for  the  virtual  potential  temperature 
and  the  speeific  humidity  show  the  effects  of  the  entrainment  very  clearly  at  the  top  of  the 
PBL  where  the  variances  are  as  much  as  50  times  larger  than  over  the  bulk  of  the  mixed 
layer.  The  larger  increase  in  the  0v  variance  over  that  of  q  near  the  top  of  the  mixed  layer 
agrees  well  with  the  results  of  Chou  et  al.  and  results  from  a  larger  normalized  virtual  poten¬ 
tial  temperature  jump  (A0v/Ov*)  than  humidity  jump  (Aq/q*).  A  comparison  of  the  LESS 
u  variance  with  that  of  Moeng’s  LES  model  shows  general  good  agreement  in  the  lower 
half  of  the  PBL  but  larger  values  for  the  LESS  results  in  the  upper  half.  The  latter  may 
be  due  to  the  secondary  maximum  in  buoyancy  at  level  4  in  the  LESS  case  resulting  in 
corresponding  increase  in  u  variance. 

Perhaps  the  most  significant  finding  in  the  research  presented  is  the  discovery  of  a 
secondary  buoyancy  maxima  at  level  4  which  appears  to  be  due  to  condensational  heating 
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which  has  resulted  trom  the  formation  of  precipitating  and  clouds  at  this  level.  While  results 
for  the  normalized  buoyancy  from  the  MASEX  and  AMTEX  studies  are  generally  in  very- 
good  agreement  with  the  LESS  results  in  the  lower  half  of  the  PBL,  neither  of  the  other 
two  studies  indicate  this  secondary  maxima  in  buoyancy  within  the  cloud  layer.  The  LESS 
results,  however,  are  in  very  good  agreement  at  least  qualitatively  with  the  model  results 
of  Deardorff.  The  vertical  profile  of  the  normalized  average  vertical  eddy  moisture  flux 
shows  good  agreement  between  the  LESS  and  MASEX  results  with  the  maximum  occurring 
for  the  LESS  at  level  2  where  both  the  w  and  q  variances  were  large. 

The  vertical  profile  of  the  normalized  average  TKE  for  the  LESS  case  study  shows 
maxima  at  levels  2  and  4.  Both  maxima  appear  to  he  primarily  due  to  the  corresponding 
maximum  in  the  horizontal  component  of  TKE  at  these  same  levels.  Both  the  horizontal 
and  vertical  components  of  TKE  are  stronger  for  the  MASEX  case  than  for  the  LESS.  This 
is  probably  due  to  the  shorter  fetch  for  the  former  resulting  in  greater  shear,  and,  hence, 
greater  turbulence.  For  both  cases,  the  horizontal  component  is  stronger  than  the  vertical 
but  just  the  opposite  is  true  of  the  vertical  flux  of  the  respective  components  of  TKE.  The 
vertical  flux  of  the  vertical  component  of  TKE  is  generally  much  stronger  than  the  vertical 
flux  of  the  horizontal  component.  For  the  LESS  and  AMTEX  data,  the  values  for  the  vertical 
flux  of  horizontal  TKE  arc  comparable,  however,  the  vertical  flux  of  vertical  TKE  is  general¬ 
ly  stronger  than  for  the  AMTEX  case  indicating  that  the  turbulent  eddies  for  the  LESS  case 
arc  more  efficient  and  capable  at  transporting  the  vertical  component  of  TKE  than  arc  the 
eddies  for  the  AMTEX.  The  analysis  of  the  ratio  of  the  vertical  flux  of  velocity  variance 
between  the  horizontal  and  vertical  components  for  all  three  case  siudies  shows  that  the  LESS 
has  a  stronger  vertical  transport  of  vertical  variance  at  the  lowest  levels  than  cither  of  the 
other  two  case  studies. 

Much  controversy  has  developed  in  the  literature  with  regards  lo  comparisons  of  verti¬ 
cal  velocity  skewness  for  observational  case  studies  and  models.  For  observational  studies 
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vertical  velocity  skewness,  Sw,  is  positive  throughout  the  Type  I  con\cetivc  PBl.  and  nearly 
constant  in  the  upper  half  of  the  PBL.  Model  Sw  is  crroneouslv  negative  in  the  'surface 
layer  for  the  surface  heating  only  simulation  probably  because  of  poor  subgrid-scale  para 
meterizations.  Unlike  the  observational  case  study  skewness,  model  skewness  continues  to 
increase  with  height  reaching  a  maximum  at  0.9  Z/Zj.  There  is  some  evidence  of  an  increase 
in  Sw  for  the  Project  LESS  at  level  4  for  the  45-60km  flight  leg.  but  this  is  probably  due 
to  the  secondary  maximum  in  buoyancy  at  this  level  discussed  earlier  in  Chapter  3. 

Comparisons  made  between  turbulence  statistics  calculated  from  the  detrended  and 
raw  data  for  TKE,  Sw  and  buoyancy  show  very  similar  profiles  for  the  first  two  quantities 
but  not  for  the  last.  Similarities  in  the  TKE  and  Sw  profiles  for  the  two  meteorological 
parameters  are  indicative  of  very  little  trend  in  the  three  wind  components.  The  difference 
in  the  buoyancy  profiles  suggests  that  a  significant  trend  at  the  synoptic  or  large-mesoscale 
level  may  exist  in  the  virtual  potential  temperature  as  evidenced  in  the  plots  of  the  raw  20Hz 
data  (Figure  2.5). 

As  for  future  work,  it  would  be  interesting  to  develop  a  prototype  Type  I  CTBL  LES 
for  which  sensitivity  studies  could  be  done  by  adding  one  physical  process  at  a  time.  Com¬ 
parisons  could  be  made  between  the  results  of  this  new  LES  study  and  those  presented  in 
this  study.  A  fellow  member  of  the  Purdue  Mesoscale  Convection  and  Climatology  (MCAC) 
research  group,  Guan-Shu  Rao,  is  actively  pursuing  such  a  project  with  Chin-Uoh  Moeng 
at  NCAR. 

Future  study  could  also  involve  comparisons  of  turbulent  statistics  bir  the  Type  I 
‘vs'  the  Type  If  CTBL  to  explore  how  the  transitions  between  the  two  CTBLs  occur  (using 
Sw  profiles,  for  example). 
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